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NOTATION 
Expanded area, zsR cdr 
h 


Ay Propeller disk area, nD*/4 

A. Fourier cosine coefficient of radial component of wake 
velocity 

A. Velocity cosine coefficient of tangential component of wake 
velocity 

A, Fourier cosine coefficient of longitudinal component of 


wake velocity 


B.. Fourier sine coefficient of radial component of wake 
velocity 
Be Fourier sine coefficient of tangential component of wake : 
velocity 4 
BY Fourier sine coefficient of longitudinal component of , 


wake velocity 


Cy Correlation allowance 

4 Elements of calibration matrix 

Con Thrust loading coefficient, T/((p/2)V,"Ap)) 

c Blade section chord length 

D Propeller diameter 

FA Froude number 

(F) nth harmonic amplitude of F 

¥5952 Force components on blade in x,y,z directions 

fu Camber of propeller blade section 

J Advance coefficient, J=V,/nD 

Jp Effective advance coefficient based on thrust identity 
Jog Effective advance coefficient based on torque identity 
Jy Ship speed advance coefficient, J=V/nD 


2.4 
Force coefficient, Fey 2/ (on D) 
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Moment coefficient, M /(on2D°) 
X,V,zZ 


} X,Y 2 
Ky Torque coefficient, Q/(pn2D°) 
Keo Centrifugal blade spindle torque coefficient, M, /(opn’D”) 
Cc 

Kp Thrust coefficient, T/(pn2D*) 
} M Moment components about x,y,z axes from loading on one 
} X,Y,zZ 

blade 
(™) nth harmonic amplitude of M 
n Propeller revolutions per unit time 
i 
E P Propeller blade section pitch 
i 
i Q Time average propeller torque arising from loading on all 
I blades, -2M 
f x 
i R Radius of propeller 
: Ro Reynolds number, cy 7V,rly 
r Radial coordinate from propeller axis 
+ Tp Radial center of hydrodynamic component of axial force, 
x 

it ae 3 ae 

4 O*H 

Radial center of hydrodynamic component of transverse 
Fy 
H force, M /(F_ ) 
H OH 

Ss Skew back of propeller blade section measured from the 

spindle axis to the midchord point of the blade section, 
| positive towards trailing edge 
1 Time average thrust of propeller, positive foward, ZF. 
t 
t Maximum thickness of propeller blade section 
t 
Vv Model speed 
f 
i Va Propeller speed of advance 
V,* Vector sum of speed of advance and rotational velocity 

at the 0.7 radius, qv? + (0. 7mnD)2)1/2 

VG, 84) Radial component of wake velocity, positive towards hub 
nth harmonic amplitude of Xe 
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> 
< 
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Tangential component of wake velocity, positive clockwise 
looking upstream for starboard propeller (left hand 
rotation), positive counterclockwise looking upstream 
for port propeller (right hand rotation) 


nth harmonic amplitude of . 
Longitudinal component of wake velocity, positive forward 
nth harmonic amplitude of Ve 


Volume mean longitudinal velocity through propeller disk 
determined from wake survey 


Taylor wake fraction determined from torque identity 
Taylor wake fraction determined from thrust identity 


Wake fraction determined from volume mean longitudinal 
velocity through propeller disk determined from a wake 
survey, (V-V)/V 


Coordinate axes 
Number of blades 


Rake of propeller blade section measured from the pro- 
peller plane to the generator line, positive aft 


Advance angle at 0.7 radius, 


ten”* [v, (r=0.07)/(0.71n »)| 


Blade strain 


Angular coordinate used to define location of blade and 
variation of loads, from vertical upward positive counter- 
clockwise looking upstream for starboard propeller (left 
hand rotation), positivé clockwise looking upstream for 
port propeller (right hand rotation), 6=~6,, 


Skew angle measured from spindle axis to projection of 
blade section midchord into propeller plane, positive 
toward trailing edge 


Angular coordinate of wake velocity, from upward vertical, 
positive clockwise looking upstream for starboard pro- 
peller (left hand rotation), positive counterclockwise 
looking upstream for port propeller (right hand rotation), 
6. .=-0 

W 


ix 


Subscripts: 
A 


Cc 


Ship to model linear scale ratio 
Kinematic viscosity of water 
Mass density of water 

Mass density of propeller blade 


Pitch angle of propeller blade section, 
1 (@/(2nr)) 


nth harmonic phase angles of F,M based on a cosine 
series (F,M)=(F,M) + spe (F,M) | cos (n0-(> 


FM 
nth harmonic Boru angles of X. based on a sine series, 


¥. cas ) + a cy. 2. sin (n6,, - (oy.*) ) 


nth harmonic Suns angles of Vv, based on a sine series, 
Lp my. )+ me Cu sin (n8,, e (oy,* * 


nth harmonic eee angles of . based on a sine series, 
x 
woe WV, i Po (V,),, sin (n0, + ($y,*) 1) 


Pitch angle of hull 


Applied values of loads 

Arising from centrifugal loading 
Value in calm water 

Arising from gravitational loading 
Arising from hydrodynamic loading 
Value of hub radius 


Indicated values of loads before calibration matrix is 
applied 


Model value 
Maximum value at any blade angular position 


Minimum value at any blade angular position 


PEAK 


SP 
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Superscripts: 


Value of nth harmonic 
Port propeller 


Peak value including variation of both time-average value 
per revolution and variation with blade angular position 


Ship value 
Value at self-propulsion point 
Starboard propeller 


Total loading from hydrodynamic, centrifugal, and 
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ABSTRACT 


Experiments are described in which the mean and unsteady 
loads were measured on a single blade of a model of the 
controllable-pitch propeller on the DD-963 Class Destroyer. 
The experiments were conducted behind a model of the DD~963 
hull under steady ahead operation, hull pitching motions, 
and simulated acceleration maneuvers. The experimental tech- 
niqves are outlined and the dynamometer and data analysis sys- 
tem described. 

The results show that ali significant loads except radial 
force are predominantly o. hydrodynamic origin. The circum- 
ferential variation of all measured components of blade load- 
ing is primarily a once-per-revolution variation, with the 
variation following approximately the variation of the tan- 
gential wake velocity. 

For sinusoidal pitching of the hull with maximum pitch 
angle of 1.85 degrees and a simulated full scale frequency 
of 0.16 hertz, the peak-to-peak circumferential variation of 
measured forces and moments increased by approximately 50 per- 
cent over the values without hull pitching. 

For simulated operation during an acceleration maneuver, 
the circumferential variation of measured forces and moments 
varied approximately as the product of ship speed and propeller 
rotational speed. At no time during the simulated acceleration 
maneuvers were the circumferential variations of loads as large 
as during full power steady ahead operation. 

For steady ahead operation, circumferential variation of 
loading determined from the model experiments agreed fairly 
well with full-scale data, but was substantially larger than 
the theoretically calculated values. 

For all conditions evaluated, the results follow close to 
previously reported results of similar experiments on a model 
of the FF-1088. 


ADMINISTRATIVE INFORMATION 
The work reported herein was funded by the Naval Sea Systems Command 
(NAVSEA 033), Task Area S0379-SLOO1, Task 19977. The work was performed 
under David W. Taylor Naval Ship Research and Development Center (DINSRDC) 
Work Unit No. 1544-296. 


The English system of units was used in the original calculations 


presented in this report. Therefore, all data are presented in the 


English units. However, the International System (SI) of metric units are 


shown in the text in parentheses following the English units. 


INTRODUCTION 


Major naval ships powered with marine gas turbines and using 


controllable-pitch (CP) propellers for thrust reversal are currently being 
added to the Fleet. Ships with CP propellers include the DD-963 Class, 
the FFG-7 Class, and the DDG-47 Class. 

Accordingly, the Navy has been conducting a research and development 
(R&D) program to establish the technology for producing reliable CP pro- 
pellers with delivered power in the range of 35,000 to 40,000 horsepower 
(26,000 to 30,000 kw) 2 As part of this program, CP propellers were 
installed on the U.S.S. PATTERSON (FF-1061) and U.S.S. BARBEY (FF-1088) 
with delivered power of 35,000 horsepower (26,100 kW). These installations 
were intended to demonstrate that CP propellers in this range of power had 
adequate reliability for application to ships with gas turbine prime 


movers. 


Because of the structural failure of the crank rings to which the 


blades of the CP propeller on the FF-1088 were bolted, R&D efforts were 
intensified. The program undertaken at DTNSRDC included: 
1. Blade Loading of CP Propellers 
a. Model measurement and theoretical prediction of blade load- 


ing on CP propellers. 


b. Model and full-scale wake measurements and theoretical pre- 


dictions of wake. 

Full-scale measurements of forces, pressures, and strains in 
CP propeller component 8. 

2. Structural Design of CP Propeller Blade Attachments. 

3. Development of Materials for CP Propeller Systems. 


The current report presents the results of work conducted under Sec- 


tion la of the CP Propeller Research and Development Program, i.e., model 


lingelo, J.J. et al, "U.S. Navy Controllable Pitch Propeller Programs," 
presented at a Joint Session of the Chesapeake Section of the Society of 
Naval Architects and Marine Engineers and the Flagship Section of the 
American Society of Naval Engineers, Bethesda, Maryland (19 April 1977). 


measurement and theoretical prediction of blade loading of CP propellers. 

Work under the other sections of this program will be reported separately. 
The present report presents experimental results obtained on a model 

of the CP propeller on the DD-963 Class Destroyer. The results of similar 


experiments on a model of the FF-1088 were reported in Reference 2. 


BACKGROUND 

Extreme care must be taken to design the blades and pitch-changing 
mechanisms of high power CP propellers so that they possess adequate 
strength including consideration of yield and fatigue stresses. This 
requires an accurate estimate of the maximum time-average and alternating 
loads under all operating conditions. High time-average and alternating 
loads occur at steady full-power ahead conditions and during high-speed 
maneuvers including full-power crash astern, full-power crash ahead, and 
full-power turns. In addition, the influence of the seaway may substan- 
tially increase the time-average and alternating loads. At present there 
appears to be no confirmed technique whereby the pertinent loads can be 
predicted to the desired accuracy. Schwanecke and Wereldsma” reviewed the 
factors affecting blade loading for propellers in general, and Rusetskiy“ 
and Hawdon et al," discussed some of the factors peculiar to blade loading 


of CP propellers. 


2eoswell, R.J. et al, "Experimental Determination of Mean and Unsteady 
Loads on a Model CP Propeller Blade for Various Simulated Modes of Ship 
Operation," The Eleventh Symposium on Naval Hydrodynamics sponsored Jointly 
by the Office of Naval Research and University College London, Mechanical 
Engineering Publications Limited, London and New York, pp 789-823, 832-834, 
(April 1976); also "Experimental Unsteady and Mean Loads on a CP Propeller 
Blade of the FF-1088 for Simulated Modes of Operation," David Taylor Naval 
Ship Research and Development Center Report 76-0125, October 1976. 


35chwanecke, H. and R. Wereldsma, "Strength of Propellers Considering 
Steady and Unsteady Shaft and Blade Forces, Stationary and Nonstationary 
Environmental Conditions," Proceedings of the Thirteenth International 
Towing Tank Conference, Report of the Propeller Committee, Appendix 2b, 
Vol. 2, pp 495-526 (1972). 


Snusetskiy, A.A., “Hydrodynamics of Controllable Pitch Propellers," 
Shipbuilding Publishing House, Leningrad (1968). 


5 


Hawdon, L. et al, "The Analysis of Controllable-Pitch Propeller Char- 
acteristics at Off-Design Conditions," Transactions of the Institute of 
Marine Engineers, Vol. 88, Series A, Part 4, pp 162-184 (1976). 


3 


Near the self-propulsion point in calm water, the time-average loads 
can probably be calculated with reasonable accuracy. However, even at 
these conditions, the variation of loads with blade angular position ap- 
parently cannot be calculated with high accuracy. Various techniques, 
including quasi-steady procedures, stripwise unstead procedures, and 
methods based on unsteady lifting surface theory, have been proposed for 
calculating the unsteady loading arising from the circumferential varia- 
tion in the inflow velocdey. However, all of these procedures require 
knowledge of the flow patterns (wake profile) in the propeller disk. In 
current practice, the wake profile is measured in the plane of the propel- 
ler behind the model hull with the propeller removed. For high-speed dis- 
placement ships of the type under consideration in this report, these 


results are usually extrapolated to full scale without making allowance for 


Sprwin, J.E., "The Development of Numerical Methods for the Computa- 
tion of Unsteady Propeller Forces," Presented to the Symposium on Hydro- 
dynamics of Ship and Off-Shore Propulsion Systems, Oslo, Norway (March 
1977). 


‘S20, T. et al, "Calculation of Unsteady Propeller Forces by Lifting 
Surface Theory," Presented to the Symposium on Hydrodynamics of Ship and 
Off-Shore Propulsion Systems, Oslo, Norway (March 1977). 


Spoddy, R.F., "A New Method for the Calculation of Unsteady Forces on 
a Marine Propeller," Presented to the Chesapeake Section of the Society 
of Naval Architects and Marine Engineers, Washington, D.C. (February 1977). 


Ivan Gent, W., “Unsteady Lifting Surface Theory for Ships Screws: Deri- 
vation and Numerical Treatment of Integral Equations," Journal of Ship 
Research, Vol. 19, No. 4, pp 243-253 (December 1975). 


10. nwanecke, H., "Comparative Calculations on Unsteady Propeller Blade 
Forces," Proceedings of the Fourteenth International Towing Tank Confer- 
ence, Report of Propeller Committee, Appendix A, Vol. 3, pp 357-397 (1975). 


a eedlin, J.P., "Propeller Excitation Theory,'' Proceedings of the 
Thirteenth International Towing Tank Conference, Report of the Propeller 
Committee, Appendix 2c, Vol. 2, pp 527-540 (1972). 


12, swell, R.J. and M.L. Miller, "Unsteady Propeller Loading - Measure- 
ment, Correlation with Theory, and Parametric Study," Naval Ship Research 
and Development Center Report 2625 (October 1968). 


(1) the change in Reynolds number and the corresponding reduction in 
relative boundary layer thickness and (2) the effect of the propeller suc- 
tion on the boundary layer and thereby the wake pattern in propeller disk. 
Although these two factors may be important for full-form ships such as 
cargo ships, they are probably not important for high-speed transom stern 
ships of the type under consideration in this report. 

Existing measurements which give information on unsteady blade load- 
ing include: 

1. Measurements of strain on the blades of the model propellers or 
full-scale propellers. However, some calculations and assumptions are 
required to convert measured strains into loads. Published data of this 
type have been summarized by Meyne.}3 

2. Measurement of bearing (shaft) forces and moments on model pro- 
pellers operating in wakes generated by model hulls or wire grid screens. 
However, these experiments give information on only some components of 
blade loading and on only those harmonics of shaft rotational speed cor- 
responding to nZ-1, nZ, and nZ+l, where n is an integer and Z is the 
number of blades. Measurements of this nature have been conducted by many 
investigators, as summarized by beésita™” and Wereldess.*” 

3. Measurements of forces and moments on individual blades of model 
propellers operating in wakes generated by model hulls or wire grid 
screens. Measurements behind model hulls have been made by Buse” and 


5 
Slauroek,~” measurements behind screens have been made by Hawdon et al, 


13 yeyne, K., "Propeller Manufacture - Propeller Materials - Propeller 
Strength," International Shipbuilding Progress, Vol. 22, No. 247, pp 77- 
102 (March 1975). 

MAvereldena, R., "Comparative Tests on Vibratory Propeller Forces," 
Proceedings of the Thirteenth International Towing Tank Conference, 

Report of the Propeller Committee, Appendix 2a, Vol. 2, pp 482-494 (1972). 
huss, E., “An Experimental Investigation of the Dynamic Forces and 
Moments on One Blade of a Ship Propeller," Proceedings of the Symposium on 
Testing Techniques in Ship Cavitation Research, The Norwegian Ship Model 
Experimental Tank, Trondheim, Norway, Publication No. 99, Vol II, pp 19- 
188 (December 1967). 


163; surock, J., "Propeller Blade Loading in Nonuniform Flow," The 
Society of Naval Architects and Marine Engineers, Propellers 75 Technical 
and Research Symposium S-4, Paper No. 4, pp 4/1-4/17 (February 1976). 


5 


and Raestad,*’ measurements in inclined flow have been made by Albrecht and 


Subcbiers-” Ketnaeeicc and Gaeeeade aad measurements on partially sub- 


merged propellers have been made by Pabay.” 


Experiments in wakes generated by screens are advantageous for 
evaluating the ability of a procedure to calculate the loading for a given 
wake since for this case, the propeller apparently does not influence the 
wake pattern. Although some good correlation has apparently been obtained 
between analytical predictions and unsteady bearing forces measured be- 


i, 2 


hind wire grid screens, correlation has been rather inconsistent 


between analytically predicted unsteady blade loads, or resulting 
strains, and measured. blade loads, or beectae os 
The mechanism by which the seaway influences the mean and unsteady 
blade loads is complex. Factors include the increased mean propeller 
loading due to increased hull resistance and the increased unsteady load- 
ing resulting from the influence of the free surface and modification of 
the flow pattern into the propeller disk. This flow pattern is influenced 
by (1) direct trochoidal velocities from the ocean waves, (2) relative 
velocities of the propeller due to ship motions, and (3) modification of 
the hull wake pattern due to the seaway and ship motions. Procedures for 


calculating the loads in a seaway are much less refined than for steady 


I psestad, A.E., "Hydrodynamic Propeller Loading in the Behind 
Condition," det Norske Veritas Research Department Report 74-31-M (1974). 
18 ibrecht, K. and K.R. Suhrbier, "Investigation of the Fluctuating 
Blade Forces of a Cavitating Propeller in Oblique Flow," International 

Shipbuilding Progress, Vol. 22, No. 248, pp 132-147 (April 1975). 


19, ednaretk, R., "Untersuchung uber die Belastungsschwankungen am 
Einzelflugel schrag angestromter Propeller," Schiffbauforschung, Vol. 8, 
No. 1/2, pp 57-80 (1969). 


0 nobay, G.F., "Time-Dependert Blade-Load Measurements on a Screw- 
Propeller,” Presented to the Sixteenth American Towing Tank Conference, 
Instituto De Pesquisas Techologicas, Marinha Do Brasil (August 1971). 


21 ereldena, R., "Last Remarks on the Comparative Model Tests on 
Vibratory Propeller Forces," Proceedings of the Fourteenth International 
Towing Tank Conference, Report of the Propeller Committee, Appendix 7, 
Vol. 3, pp 421-426 (1975). 
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operation in calm water. Tapaki7” gives a good review of the mechanisms 
and procedures for predicting the effect of the seaway on bearing forces 
which, in principle, also apply to unsteady loading on an individual 
blade. Keil et a4? and Watanabe et ai.** present strain measurements 
on the blades of full-scale propellers in both calm and rough seas. 

Apparently no rational analytical procedures are available for 
accurately calculating the time-average loads per revolution or the un- 
steady loads including variation with blade angular position during crash- 
ahead or crash-astern maneuvers. These loads may depend on many factors 
including the time rate of change of propeller pitch p (for CP propel- 
lers), time rate of change of rotational speed fh, time rate of change of 
ship speed V, propeller blade-section stall, cavitation, ventilation, 
flow separation from the hull, and large interactions between the propel- 
ler and the hull. Some of these factors are discussed and considered by 
Hawdon et a 

For turns, the factors affecting the time-average loads per revolu- 
tion and the unsteady loads are somewhat the same as those affecting the 
loads under crash-ahead and crash-astern conditions except that for turns, 
there is a relatively large drift angle of the flow into the propeller. 
This drift angle tends to increase the circumferential nonuniformity of 
the flow into the propeller, thereby increasing the unsteady loading. 
However, this circumferential nonuniformity of the inflow tends to be off- 
set by the lower values of ship speed and propeller rotational speed in 


turns compared to steady ahead operation. 


Prior to the present R&D investigation, no experimental measurements 


existed to the authors' knowledge which showed the time-average loads and 


a2 agaki, R., “Propulsion Factors and Fluctuating Propeller Loads in 


Waves," Proceedings of the Fourteenth International Towing Tank Confer- 
ence, Report of Seakeeping Committee, Appendix 7, Vol. 4, pp 224-236 
(1975). 


2S yeti, H.G. et al, "Stresses in the Blades of a Cargo Ship Propeller," 
Journal of Hydronautics, Vol. 6, No. 1, pp 2-7 (January 1972). 


24 uatanabe, K. et al, "Propeller Stress Measurements on the Container 
Ship HAKONE MARU," Shipbuilding Research Association of Japan, Vol. 3, No. 
3, pp 41-51 (1973). 
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circumferential variation of loads with blade angular position on CP pro- 
pellers behind a hull under a wide range of operating conditions. An ex- 
perimental program was therefore undertaken to measure the six components 
of loading (Figure 1)* on model CP propellers operating behind model hulls. 
The initial experiments were conducted on a single-screw ship, namely, a 
model of the FF-1088. These results were reported previously." The sec- 
ond set of experiments were conducted on a twin-screw ship, namely, a 
model of the DD-963 Class. These results are presented in the current 
report. 

For the DD-963 Class, the experimental conditions included 
(1) steady-ahead operation near the self-propulsion point, (2) steady- 
ahead operation near the self-propulsion point with forced dynamic pitch- 
ing of the model hull, and (3) simulated acceleration operation. 

Results for the steady ahead operation were correlated with predic- 
tions based on unsteady lifting surface theory as developed by Tsakonas 
et Oi and with the quasi-steady method of Mccarthy ,7° and with strains 
measured on the full-scale propeller. 

Blade loading measurements were made on the propeller on the star- 
board shaft since this shaft had a larger rake angle than the port shaft. 
The propellers used in these experiments were DINSRDC propellers 4660 
(right hand rotation on port shaft) and 4661 (left hand rotation on star- 
board shaft), which were made of aluminum; see Figure 2 and Table 1.** 
The hull of the DD-963 Class was represented by DINSRDC model hull 
5265-1B; see Figure 3. 


*Figures are presented following the section on acknowledgments. 


! *kThe tables are presented following the figures. 

| 

: *Steakonas, S. et al, "An Exact Linear Lifting Surface Theory for 
Marine Propeller in a Nonuniform Flow Field," Journal of Ship Research, 
Vol. 17, No. 4, pp 196-207 (December 1974). 


26uccarthy, J.H., "On the Calculation of Thrust and Torque Fluctuations 


of Propellers in Nonuniform Wake Flow,'' David Taylor Model Basin Report 
1533 (October 1961). 


EXPERIMENTAL TECHNIQUE 


FACILITY AND DYNAMOMETRY 

All experiments were conducted on DTNSRDC Carriage I using basically 
the same dynamometry and hardware as previously described in Reference 2. 

The port propeller, on which blade loads were not measured was driv- 
en from inside the model hull as would be the case in a self propulsion 
experiment. The propeller rotational speed, which could be controlled 
independently of the starboard propeller, was measured via a toothed- 
pickup and recorded on a digital voltmeter. The time-averaged thrust and 
torque were measured for selected runs by a transmission dynamometer. 

The starboard propeller, on which blade loads were measured, was 
located in its proper position relative to the model hull but was isolated 
from the hull and driven from downstream (see Figure 4). This downstream 
drive system was necessary in order to obtain the required characteristics 
of the system for measuring unsteady loading. The general criteria for 
the design of an unsteady force measuring system are: 

1. The support structure of the force measuring system should be 
soft mounted and possess a large mass to eliminate transmission of extra- 
neous vibration to the system. 

2. The natural frequency of the system should be well above the 
highest frequency of the quantities to be measured (to avoid phase shift 
and amplification of the signal). 

3. The system response in the force magnitude range should be suf- 
ficiently large to be measurable (sensitivity). 

4. The system should be free of interaction, that is, each measur- 
ing element should respond only to that force or moment which it is in- 
tended to measure. 

These four major aims are not complementary. The high natural fre- 
quency requires a stiff, rigid system whereas high sensitivity requires 
an elastic, soft system. The necessary compromise results in some inter- 
action between the force-measuring elements. 

Criterion 1 dictated that a massive flywheel be used, and Criterion 
2 dictated that this flywheel be connected to the sensing elements (locat- 


ed inside the propeller hub) by a short thick shaft. Therefore, because 


of the geometry of the hull and shafting of the configuration under eval- 
uation, it was not feasible to achieve both these criteria with an up- 
stream drive system from inside the model hull. Criteria 1 and 2 con- 
trolled the minimum allowable beam and draft of the downstream body and 
the maximum allowable clearance from the bow of the downstream body to the 
propeller. Although the downstream body may exert some influence on the 
flow into the propeller, that location was considered necessary in order 
to meet these measuring criteria. The influence of the downstream body 

on the flow into the propeller is discussed in the section on experimental 
results. 

The drive and mounting system was basically the same as that used 
in the DTNSRDC BASS dynamometer which has been described by irendaa. 
Utilized from this dynamometer were the propeller (tail) shaft, drive 
shaft with flywheel, belt-type (quiet) transmission, and sliprings. 

Power to rotate the propeller was supplied by a d-c permanent-magnet 
servomotor capable of delivering up to 33 foot-pounds (45 N-m) of torque. 
The electrical power to this motor was delivered through a precision 
solid-state motor controller so that the shaft revolution rate could be 
controlled very accurately and held over the wide range of propeller 
torque loadings required for some of the experimental conditions. Mounted 
on the propeller shaft was a digital encoder that generated electrical 
pulses as a function of shaft angular position. Two types of pulses were 
generated: a single pulse per revolution and a multipulse per revolution 
(90 equally spaced pulses for the current experiment). The single pulse 
was syncronized with the reference line of the instrumented propeller 
blade. The pulses generated by this encoder are accurate to within 0.01 
degree. 

The downstream body which housed the drive system was basically that 
used by Dobay”” but modified to allow deeper submergence and an inclined 
shaft angle. Both the body housing (the drive system was soft mounted to 
this body) and the model hull were rigidly attached to a pitch-heave 


27) seandau, J.H., "Static and Dynamic Calibration of Propeller Model 
Fluctuating Force Balances," David Taylor Model Basin Report 2350 (March 
1967); also Technologia Naval, Vol. 1, pp 48-74 (January 1968). 


10 


oscillator which, in turn, was rigidly mounted on the towing carriage. 
This arrangement enabled che model hull and the drive system to be dy- 
namically pitched together while maintaining independent support from one 
another. 

The sensing elements were flexures to which were bonded high- 
sensitivity, semiconductor strain-gage bridges. The design of these 
flexures has been described by Pobey.-” There were three flexures, each 
of which measured two components of blade loading. Flexure 1 measured 
components F. and M3 Flexure 2 measured components es and M, and Flexure 
3 measured components F and M, (Figures 1, 5, and 6). An arrangement of 
three separate flexures rather than one to measure all components of blade 
loading was adopted because it appeared to result in higher natural fre- 
quencies (Criterion 1), higher sensitivities (Criterion 3), and lower 
interactions (Criterion 4), than would have resulted had a single flexure 
been used. 

The flexures were mounted inside a propeller hub which was specifi- 
cally designed for these experiments (Figure 6). Only one flexure could 
be mounted at a time, because of space limitations, and this necessitated 
duplicate runs, as discussed later in the section on experimental condi- 
tions and procedures, 

The strain-gage bridges were excited by a common d.c. voltage 
source, transmitted through the sliprings on the propeller shaft. The 
constant-current excitation used by Dobay”” was not employed in the pre- 
sent experiment because it appeared to be too sensitive to temperature. 

The voltage output from the flexures (due to blade loading) was 
transmitted through the sliprings to individual amplifiers (NEFF 119-121). 
These amplifiers utilized field effect transistors to produce an extremely 
high input-impedance (100 megohms, minimum). This high impedance essen- 
tially eliminated slipring noise to the amplifier. The voltage signals 
were transferred across the sliprings in the presence of only a small 
amount of noise-producing current. The amplifiers used here had zero- 
phase shift qualities in the d.c. to 20 kilohertz range. They were chopper- 
stabilized to enable both the steady and unsteady signals to be recorded 
simultaneously. This signal-conditioning system was essentially the same 


as that used by Dobay. 7? 


11 


al ee ee eee eal 


—7= 


The signals were then digitized and analyzed by using a Model 70 
Interdata Digital Computer, and then stored in digital form on a nine- 
track magnetic tape. The on-line analysis of the data is discussed in the 


section on data acquisition and analysis. 


CALIBRATION 

Prior to the experiment, each flexure was statically calibrated in 
air to establish flexure sensitivities, interactions, and linearity over 
the loading range of interest. These calibrations were conducted with the 
flexures mounted in a calibration stand, with the flexure electrical ca- 
bles connected through the flywheel and drive assembly as in the experi- 
ment. Each flexure was subjected to independently controlled forces in 
the axial, transverse, and radial directions (i.e., a ee and Fa re- 
spectively) and to independently controlled moments about the axial, 
transverse, and radial directions (i.e., M,» Ms and M,> respectively); 
see Figure l. 

The static calibration showed that all flexures had a linear re- 
sponse over the load range of interest. Table 2 shows the interaction 
matrix. These calibrations indicated that all flexures had good sensi- 
tivity except Re whose sensitivity was slightly lower than desirable. 

The interactions were small except for the effect of M, on Fo The infe- 
rior characteristics of the rs flexure is not considered a serious short- 
coming since ¥ arises primarily from centrifugal loading and can be 
analytically calculated. In addition, no significant variation of Fo 
with blade angular position was anticipated. Flexure 3, which measured 
‘. and Mo» was further evaluated by correlation of air-spin experiments 
with analytically calculated results, as discussed later. The interac- 
tions were taken into consideration during data analysis. 

The flexures used in this experiment had been dynamically cali- 
brated by Dobay”? to determine the frequency range over which unsteady 
forces and moments could be reliably measured. In this procedure, an 
electromagnetic shaker in air was used to apply a relatively constant, 
maximum amplitude, variable-frequency force or moment-excitation in all 
six-component directions to all six flexure elements. The force or mo- 


ment amplitude imposed by the shaker was monitored through an extremely 
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light-weight, strain-gaged single flexure element. The measured lowest: 


natural frequencies of the three flexures in air were as follows: 


Frequency (hertz) Mode 
Flexure 1 550 


Flexure 2 450 
Flexure 3 282 


The measured amplification factor (ratio of output amplitude to in- 


put amplitude) and phase shift for all three flexures was as follows: 


Frequency Range (hertz) 0 to 60 60 to 120 
Phase Shift (degrees) 0 to 0.05 0.05 to 0.15 
Amplification Factor 1.00 1.00 to 1.05 


In the previous Sxperinent the natural frequency of the flexures 


was found to be reduced by 55 percent when submerged with blades attached. 
As a result, it was concluded that the flexures had a "true" dynamic re- 
sponse up to at least the third harmonic and no greater than a five per- 
cent amplification up to the sixth harmonic. Because the blades used on 
the present experiment were lighter and smaller than those used on the 
previous experiment, it was assumed that the natural frequency of the 
flexures would be reduced to a lesser extent when submerged with blades, 
so that the dynamic amplification would be less. This assumption was 
supported by the increase in frequency of the extraneous signals appearing 
in the unfiltered experimental data as discussed in the experimental re- 
sults section. 

The propeller shaft drive and soft-mount support system were dynam- 
namically loaded in the vertical, longitudinal, and transverse directions 
to obtain the lowest natural frequencies of the system. The natural fre- 


quencies of the system in air were found to be: 


Mode Natural Frequency (hertz) 
Vertical bending nae 


Horizontal bending 6.0 
Axial 4.6 


The support system had a low resonant range; however, the soft- 


Mount system was specifically designed to prevent towing-carriage 
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oscillation (with the resonance at 100 to 200 hertz) from being trans- 
mitted to the blade flexures. Based on the measured resonance, it is 
concluded that the soft-mount system should successfully meet this ob- 
jective. Although some resonances were close to the propeller rotational 
speed for some experimental conditions, it was considered more desirable 


to isolate the system from towing-carriage vibration. 


EXPERIMENTAL CONDITIONS AND PROCEDURES 

Experiments were conducted at several conditions including steady- 
ahead operation, simulated pitching of the hull, and simulated accelera- 
tion. All conditions were run with the model hull rigidly attached to its 
support, with no freedom to sink or trim, and with essentially equal rota- 
tional speed on the port and starboard propellers. 

The steady-ahead condition is defined in Tables 3 and 4. The simu- 
lated full scale ship speed and propeller rotational speed for this condi- 
tion were determined from model self-propulsion data* at simulated dis- 
placement of 7,800 tons (7,920 tonnes) including corrections for wind drag 
at zero true wind and a three-percent margin on effective power with 
Com 0.0005. 


The trim and draft at this speed were obtained from Reference 28. 


These had been determined by setting the specified still water trim (even 
keel) and draft (19.5 feet (6.40 m) full-scale equivalent), attaching 
the model to the carriage so that it was free to trim and sink, running at 
the specified speed, and locking the model at this equilibrium trim and 
draft. The equilibrium sinkage was 0.5 feet (15.4 cm) at the bow and 3.0 
feet (98.4 cm) at the stern. 

Runs simulating hull pitching were conducted at the same conditions 
as the steady-ahead run, except that the hull pitch was varied. Two types 
of runs were conducted: (1) quasi-steady simulation in which the hull 


pitch angle was set at various fixed positions, and (2) unsteady 


*DINSRDC experiments 21 and 22 on Model 5265-1B. 


2855, W.G., "The Effect of Speed on the Wake in Way of the Propeller 
Plane for the DD-963 Class Destroyer Represented by Model 5265-1B,"" David 
Taylor Naval Ship Research and Development Center Report SPD-311-37 (July 
1975). 
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simulation in which ) was varied sinusoidally with time. For the quasi- 
steady simulation, runs were conducted at five different values of », from 
1.85 degrees bow up from the calm water equilibrium ¥W=V oy) to 1.85 de- 
grees, bow down from Yow (Tables 3 and 4). For the unsteady pitch simu- 
lation, the value of y was varied sinusoidally about Yow with an ampli- 
tude of 1.85 degrees and a frequency of 0.8 hertz.* The selected scaled 
amplitude and frequency were within the predicted response characteristics 
of the DD-963. All runs were conducted in calm water; therefore, the re- 
sponse of the hull to the seaway was simulated but the seaway was not 
simulated. 

Acceleration runs were conducted based on analytical dynamic simu- 
lation studies of the DD-963. 7? The experimental conditions followed run 
7501061 of Reference 29, which was an acceleration from 8.7 knots to full 
power (see Tables 3 and 4 and Figure 7). Trim and displacement were fixed 
at the values corresponding to the self-propulsion condition (Condition 1 
of Table 3). Two types of runs were conducted: (1) quasi-steady runs in 
which all quantities including model speed V, rotational speed n, and pro- 
peller pitch P were held constant (V=n=P=0), and (2) unsteady runs in 
which V was varied with time but n and P were held constant (V>0, n=P=0). 
For the quasi-steady simulation, runs were conducted at five different 
combinations of V, n, and P. The conditions for each run represent the 
conditions at one instant of time during a "true" acceleration in which 
V, n, and P vary with time. Thus, one "true" acceleration run is repre- 
sented by five steady runs which do not simulate the time rate of change 
of V, n, and P. For the unsteady simulation, runs were conducted at the 
same five combinations of fixed n and P as used for the quasi-steady simu- 
lation, and V was varied with time (the same variation was used for each 
run) representing an acceleration of the model hull (Figure 7). For each 
of these runs, data are of interest only near that value of V which oc- 


curred concurrently with the fixed values of n and P during the "true" 


* 
Full scale equivalent frequency is 0.16 hertz. 


2 nobis, C.J. and T.R. Harper, "Propulsion Dynamics Simulation of the 
DD-963 Class Destroyer," Propulsion Dynamics, Inc., Report 74R1B (January 
1975). 
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acceleration (V>0, n#0, P#0). Thus, one "true" acceleration run is rep- 
resented by five runs which simulate the proper time rate of change of V 
7 but not the proper time rate of change of'n and P. The quasi-steady and 
unsteady acceleration simulations were for the same condieiege, the only 


difference being that V=0 for the quasi-steady simulation whereas V>0 for 


the unsteady simulation. In general, P varies with time during a "true" 
acceleration run; however, for the acceleration run under simulation here, 
P was constant throughout the portion of the run simulated. 

For the unsteady acceleration runs, the carriage speed was manually 
varied with time in a carefully controlled manner. This was achieved with 
the aid of an inked pen on a two-dimensional Cartesian plotter. In one 
direction, the pen was controlled so that it moved linearly with time, 
and in the orthogonal direction, it was controlled so that it varied with 
the instantaneous carriage speed. When an acceleration maneuver was to be 
executed, the switch moving the pen with time was turned on and the car- 


riage operator manually varied the carriage speed so that the inked pen 
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followed a prescribed velocity versus time curve. 

As discussed earlier, each of the three load-sensing flexures 
measured only two components of blade loading. Therefore, each of the 
experimental conditions described in Table 3 was run with each of the 


three blade loading flexures. 


The blade pitch was set by using a Sheffield Cordax 300 measuring 
machine. In order to change either the blade pitch or the flexure, the 
propeller had to be removed from the drive system. 

Air-spin experiments were conducted with all three flexures over a 
range of rotational speeds in order to isolate the effects of centrifugal 


and gravitational loading from hydrodynamic loading. Supplemental experi- 


: 
: 
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ments were conducted to assess the influence of the downstream dynamo- 


meter boat on the flow in the propeller plane. These supplemental exper- 
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iments consisted of wake surveys in the propeller plane at the self- 
propulsion point (Condition 1 in Table 3) with and without the downstream 
body, but without the propeller. These wake surveys yielded a direct 


measure of the change in the velocity distribution through the propeller 


disk attributable to the downstream body. The details of these wake sur- 
veys will be reported in a future DTNSRDC report.* 


DATA ACQUISITION AND ANALYSIS 

Data were collected, stored, and analyzed on-line by using a Model 
70 Interdata Digital Computer. A special-purpose computer program was 
written with options for analyzing each of the three basic types of runs: 
(1) steady ahead, (2) dynamic hull pitching, and (3) unsteady acceleration. 
These types of runs have already been discussed in detail. 

The program allowed the propeller blade force and moment data to be 
sampled and stored on magnetic tape as a function of shaft position. 
Sampling was triggered by external pulses generated by a Baldwin Digital 
Encoder mounted on the propeller shaft, as discussed earlier. Pulses were 
generated as a function of shaft angular position; hence, the sampling of 
blade force and moment data was related to shaft position. There were two 
cdtpute from the shaft encoder; a single pulse per revolution and multi- 
pulse (90 pulses per revolution for the current experiments). 

When the experimental condition was achieved, the computer operator 
initiated the data collection cycle. The program "waited" until the sin- 
gle pulse occurred; when the single pulse occurred, the computer again 
"waited" for the occurrence of the first following pulse of the 90 pulses; 
date were then sampled for all channels through an analog-to-digital con- 
verter and stored in computer memory. This process was repeated for 180 
pulses, or two shaft revolutions. At the same time, the program "read" 
two frequency counters into core memory which measured model velocity V 
and propeller rotational speed n. The values of V and n were measured by 
counting the pulses from geared wheels attached to the towing carriage 
drive system and to the propeller shaft, respectively. The values of V and 
n were averaged over two shaft revolutions. Thus, there was an average V 
and n corresponding to each pair of two consecutive revolutions. 

After two revolutions of data were sampled and stored in core memory, 


the data were transmitted from core to a nine-track digital tape recorder. 


i ont 


The transfer time was small and no pulses were missed during the transfer. 


* 
These wake surveys were conducted by R.F. Roddy, DINSRDC Code 1524. 
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The data collection cycle proceeded continuously until the operator dis- 
engaged the computer. The sampling procedure was the same for all t: pes 
of experimental conditions, and at the completion of an experimental run, | 
all data were stored on magnetic tape and were available for analysis im- 
mediately or at any later time. For the analysis, the computer operator 
selected the appropriate option of the program depending on the type of 
run, i.e., (1) steady ahead, (2) dynamic hull pitching, or (3) unsteady 
acceleration. 
The appropriate calibration factors were stored in the computer and 
considered in the analysis. However, since only two of the six components 
of blade loading were measured during a given run, the interactions between 
the various loading components could not be considered during the on-line 
analysis. The interactions were taken into account later after measure- 
ments were completed with all three flexures for a given condition. 4 
For the steady ahead condition, blade force and moment data at each : 
4-degree increment of blade angular position were averaged over the number 
of cycles recorded (usually over more than 200 cycles). Spurious data not 
related to shaft position are averaged out by this method. An harmonic 
analysis was then performed on the average wave forms of the blade loading 
components. This gave the amplitude and phase of the first 16 harmonics. 
For the dynamic pitch runs, the hull pitch angle varied sinusoidally 
with a frequency of 0.8 hertz. A position potentiometer translated bow 
vertical displacement into hull pitch angle, and this was read into the 
computer in the same manner as blade loading components. During dynamic 
pitching, the shaft rotated independent of the pitch oscillator. During 
a single propeller revolution, 90 pitch positions were measured. Thus, 
to correlate pitch angle position and revolution, an average pitch must be 
taken over each revolution. 
Fourteen dynamic pitch angle positions were selected for analysis. 
These were characterized by pitch angle y and the sign of the time rate of 
change of pitch angle ». The computer calculated the average yj and sign 
of } corresponding to each propeller revolution. Based on these calcu- 
lated average values of y and sign of ~, each propeller revolution was 
either placed in a suitable hull-pitch angle category or discarded if its 
average ~ fell outside the tolerance band of all the 14 specified values 
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of ~. Several passes down the towing tank were required in order to obtain 
a sufficient number of samples. After all the data had been sorted based 
on w, sign of }, and tolerance, the cycles for each combination of , and 
sign of } were analyzed in exactly the same manner as the data for the 
steady-ahead condition at fixed yp. 

For unsteady-acceleration runs, the model speed V varied with time 
t. During an acceleration run, data, including a measure of V, were sam- 
pled and stored in the same manner as for the steady-ahead runs. 

Five values of V were specified for analysis. For each acceleration 
run and for each specified V, the computer selected the propeller revolu- 
tion which had the average value of measured V nearest to the specified V. 
However, because only one revolution at each specified velocity was ob- 
tained for a single acceleration run, each such run was repeated five 
times. This yielded five revolutions at each specified velocity. All the 
cycles for each specified V were then analyzed in exactly the same manner 
as the data for the steady-ahead conditions. 

Thus, the on-line analysis system yielded average wave forms and 
harmonic analyses of the average wave forms for steady-ahead conditions, 
for specified conditions of ~, sign of ~ during the dynamic pitch cycle, 


and for specified velocities V during the acceleration operation. However, 


these on-line results are preliminary because: 

1. They do not consider the interactions between the various load 
components. These interactions were determined during the static cali- 
bration of the flexures. 

2. They include the complete measured signals with no filtering. 
As discussed in the section on experimental results, some extraneous sig- 
nals near the natural frequency of the flexure being used appeared to be 
superimposed on the signals generated by blade loading. 

3. They include the effect of centrifugal and gravitational load- 
ing on the aluminum model propeller. The centrifugal and gravitational 
components of loading were measured separately during air-spin experiments, 
as discussed earlier. 

4. They do not have any corrections for the influence of the dyna- 


mometer boat. These corrections are discussed later. 
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5. The bending moments were calculated about the radial location of 
the strain gages on the flexures, rather than about the shaft axis or some 
desired radius on the blade. 

Final analyses were conducted after completion of the experiment to 
consider interactions, to filter out extraneous high frequency noise, to 
isolate hydrodynamic loading from centrifugal and gravitational loading, 
to correct for the dynamometer boat, and to resolve bending moments as 
desired. These analyses were conducted using a Control Data Corporation 
(CDC) 6700 Computer. 

For each condition, the average wave form for each of the six load- 
ing components was multiplied by the inverse of the calibration matrix 


given in Table 2. 
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This matrix multiplication was performed at 4-degree increments of blade 
angular position. 

An harmonic analysis was then performed on the signals corrected for 
the interactions. Based on an harmonic analysis of the wake in the pro- 
peller plane, °° it was judged that there should be no significant loading 
of hydrodynamic origin at frequencies above ten times shaft frequency. 
Therefore, the wave form was then reconstructed by using the first ten 
harmonics of shaft frequency. This reconstruction using only the first 
ten harmonics had the same effect as filtering out all frequencies above 


ten times shaft frequency. 


Wouemings , D.E., "Numerical Prediction of Propeller Characteristics," 
Journal of Ship Research, Vol. 17, No. 1, pp 12-18 (March 1973) 
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Corrections were made to the mean values of the measured loading 
components to account for centrifugal loads and the influence of the dy- 
namometer boat, and to the first harmonic of the measured loading compo- 
nents to account for gravitational loads. The derivation of these correc- 
tions is discussed later. 

From the known values of the three measured force components and 
three measured moment components, the values of the bending moments about 
the shaft centerline and bending moments normal to the nose-tail line at 
the 0.3 and 0.4 radii were calculated. These bending moments were calcu- 
lated at every 4 degrees of blade angular position, and harmonically ana- 
ly’ :d. The wave form was reconstructed by using the first 10 harmonics of 
blade angular position, in exactly the same manner as was used for the 
other components of blade loading. 

Plots of the data were generated by the CDC computer system using a 
Calcomp Plotter. 


ACCURACY 


During the experiments for steady operation, V=0, and dynamic 


pitching, #0, where many revolutions of data were averaged during a sin- 
gle run, the standard deviations of speed V, rotational speed n, forces, 
and moments were computed, assuming the distribution in these variables at 
a given condition follows the normal probability distribution. For forces 
and moments, the standard deviation was calculated at every increment of 
blade angular position at which forces and moments were recorded. An 
error band around the data mean was then represented using the standard 
deviation multiplied by a factor dependent on the confidence level chosen. 
For the present analysis, the factor of 1.96 was selected which corre- 
sponds to a confidendence level of 95 percent. A confidence level of 95 
percent indicates a confidence (or probability) that 95 percent of the 
data considered falls within the error band. For a given run the average 
error (95 percent confidence band) in model speed V was approximately 
+0.005 foot per second (1.6 mm/s), while the error in rotational speed n 
was less than 0.001 revolution per second. The very low error in n re- 
sulted from the use of a precision solid-state motor controller as dis- 


cussed in the section on facility and dynamometry. 
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For a given steady run (V=0, ~#0) the error (95 percent confidence 
band) in measured forces and moments, fluctuated from +5 to +10 percent of 
the mean signal, depending on the circumferential blade position. Figure 
8 demonstrates the variation in error in one revolution of the uncorrected, 
raw Fy signal at the full power condition. This figure represents the 
general trend of all the force and moment components measured. 

Besides the fluctuation in signals occurring in a given run, the 
overall accuracy of the data is dependent on the repeatability from one 
run to the next. An effort was made to set experimental conditions iden- 
tically on repeat runs; however, the propeller rotational speed and model 
velocity were set by hand, so some variation was unavoidable. Table 5 
demonstrates the variation in the measured experimental conditions and the 
raw data for the FY component for 11 repeat runs. The variation in the 
mean force was +4 percent over all the runs, but on a given day the vari- 
ation averaged +2 percent. The same trend can be observed in rotational 
speed, model velocity and the harmonic force components. This day-to-day 
variation could be due to different operators setting the experimental 
conditions, slight variations in the draft of the model, and variations in 
the gain of the sensing electronics. The variations shown for t are 
typical of all the measured force and moment components. 

For all experimental conditions the rotational speed of the port and 
starboard propellers were intended to be equal. However, some exploratory 
runs were conducted to determine whether the mean or unsteady loads, which 
were measured on the starboard propeller, were influenced by the rotational 
speed of the port propeller. At a fixed value of rotational speed on the 
starboard propeller ny» the rotational speed on the port propeller n_ was 
varied. The results showed that there was no measurable effect of n_ on 
the mean or unsteady loads in the region 0.95 n. a 85 & 1.05 n,- For all 
runs for which data are presented, 0.99 n. = Bs 1.01 nes therefore, the 
results presented are not measurably influenced by inaccuracies in = 

For the unsteady acceleration (0), the average of the five values 
of V and n for which data are presented during the unsteady runs was 
generally within +0.2 foot per second (6.5 cm/s) and +0.2 revolution 


per second of the target values respectively. 
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For runs with fixed hull pitch angle y, (=0), the value of wy could 
be controlled to within +0.005 degree. For dynamic pitch runs +0, the 
selection of a propeller revolution at a specified » necessitated a toler- 
ance of 0.1 degree to ~; however, the average value of y for which data 
are presented during the unsteady runs was generally within 0.02 degree of 
the target y. 

Considering all sources of error including deviations during a run 
and inaccuracies in setting conditions, the model scale forces and moments 
presented in this report are generally considered to be accurate to within 


(plus or minus) the following variations: 
M 


in-lb (N-m) 


Steady ahead V=0,)=0 | 0.1 (0.4) | 0.2 (0.9)| 0.2 (0.02)| 0.4 (0.06) 
Dynamic pitch V=0,¥#0] 0.2 (0.9) 10.3 (1.3)| 0.4 (0.04)| 0.6 (0.07) 
Acceleration V>0,¥=0 | 0.3 (1.3) ]0.4 (1.8)] 0.6 (0.07)}] 0.8 (0.09) 


The values are somewhat more accurate for the steady-ahead runs than 
for the time-dependent runs, because the experimental conditions could be 
controlled more precisely for the steady runs and the measured forces and 
Moments were averaged over many more revolutions of the propelier. The 
time-average values per revolution (based on 90 samples per revelolution) 
are slightly more accurate than the maximum values (based on one sample 
per revolution) which took into account the variation with blade angular 
position. Further, the peak values may have been slightly influenced by 
the dynamic response of the flexures, as discussed in the section on cali- 


bration. 
EXPERIMENTAL RESULTS 


LOADING COMPONENTS 

The basic loading components are shown in Figure 1. For a right- 
hand propeller the sign convention follows the conventional right-hand 
rule with a right-hand Cartesian coordinate system. For a left-hand pro- 
peller all the loads are the same, but for this case the sign convention 


follows a left-hand rule with a left-hand Cartesian coordinate system. 
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The sign conventions for both right-hand and left-hand propellers are 
shown in Figure 1. In all pertinent figures and tables throughout this 
report the blade loading components are listed in the following order: 
1. Components measured by Flexure 1: 
a. Fy - axial force, or thrust per blade. 
b. * - bending moment about the axis normal to the shaft axis 
f at r=0. This moment is generated primarily by the eS component of 
force. 
2. Components measured by Flexure 2: 
a. Fy - tangential force, or force normal to the propeller axis 
and the spindle axis. 


b. M, - moment about the propeller axis, (r=0), or torque per 


0 eewepemesce pir See oe ee ee 


blade. This moment is generated primarily by the Fy component of 
| force. 
3. Components measured by Flexure 3: 

a. F. - radial force, or force parallel to the blade spindle 
axis. 

b. M, - moment about the spindle axis, or spindle torque. 
4. Supplemental components which were derived from the components 
listed above (not derived for all conditions). 

a. Mo. 3 = FY. oo -0.3R) cos $9.3 + Fy a -0.3R) sin $9 3 


- bending moment applied on the spindle axis about the axis inter- 
secting the spindle axis at r=0.3R and parallel to the expanded 


pitch line at r=0.3R. The My 3 vector as defined by the convention- 


al right-hand rule for a right-hand propeller (left-hand rule for 
left-hand propeller) intersects the plane normal to the propeller 
axis at the angle > 37 tan” "(Py 3/0-3 ™D) and is directed so that 


a positive M puts the face (pressure side) of the blade in 


0.3 
tension. 


b. a iin. -0.4R) cos 0.4 + ay vieatene -0.4R) sin 9.4 


- bending moment applied on the spindle axis about the axis inter- 


By * 


secting the spindle axis at r=0.4R and parallel to the expanded 
pitch line at r=0.4R. The My 4 
tional right-hand rule for a right-hand propeller (left-hand rule 


vector as defined by the conven- 
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for the left-hand propeller) intersects the plane normal to the pro- 
peller axis at the angle $6 “ce tan” (P, z/(0.4mD)) and is directed 
so that a positive My 4 puts the face (pressure side) of the blade 


in tension. In calculating M from the experimental values of the 


0.4 
three measured forces and three measured moments, an adjustment was 
necessary to allow for the contribution of loading in the region 


0.4R>r>r, =0. 3R where r, is the hub radius. It was estimated that for 


all harmonics fetes, the time average values, 3 percent of MY and 

mM, was contributed by the loading in the region 0.4R>r>r, . _— 

estimates were based on a refinement of the method of Cummings for 

the time-average values, and the method of Tsakonas et ai*? for the 
unsteady values. 

Hydrodynamic, centrifugal, and gravitational loads may contribute to 
each of these components of loading; however, for some components the 
centrifugal loads and/or gravitational loads are zero, as discussed in 
the section on centrifugal and gravitational loads. 

Each component of loading is generally presented as a variation of 


the instantaneous value with blade angular position 6 and as a Fourier 


series in blade angular position in the following form: 


F.M(6) = (FM) + 20, (FM), cos (n® - (¢p aq) 


where F,M = circumferential average value of F,M 
(F,M) = amplitude of the nth harmonic of F,M 
@ = angular position about the propeller axis, positive counter- 
clockwise from the vertical upward looking upstream for 
starboard propeller (left-hand rotation) positive clockwise 
looking upstream for port propeller (right-hand rotation) 


(op wn = phase angle of nth harmonic of F,M 

9 

where the reference line of the tlade is the spindle axis; see Figure 2 
and Table 1. 


The components My and My are the most important for determination of 


the time-average and unsteady stresses in the hub mechanism of an actual 


Or ST RT 
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controllable pitch propeller. The components M and M are the most 


0.3 0.4 
important for determination of the time-average and unsteady stresses in 


the blades of a propeller. 


CENTRIFUGAL AND GRAVITATIONAL LOADS 

The results of the air-spin experiments, corrected for interactions, 
are presented in Table 6. The time-average values arise from centrifugal 
force whereas the first harmonic arises from gravitational force. There- 
fore, the mean values should vary as x where n is the propeller rotation- 
al speed, and the first harmonic should be independent of n. 

For the mean values, which arise from centrifugal force, signficant 
nonzero values were obtained only for the ta Fr 7 and M, components. 
Any realistic propeller would have nonzero values of centrifugal loading 
components me and Ri Nonzero values of Le and se are pro- 
duced by the nonzero values of skew and rake of the propeller evaluated. 
The components 1 and ). should be zero for any geometry, however, 

a small value of (Fe was measured. This small nonzero Fo). probably 
arises from inaccuracies in the air-spin experiment and interaction 
matrix. For all components the experimentally determined mean value 
varies essentially as a. The experimental air-spin results were faired 
so that the values of the mean loading components used for separating 
hydrodynamic loads from total loads varied exactly as a” 

For the first harmonic loads, which arise from acceleration due to 
gravity, nonzero values were obtained only for the i M,. and Fo flexures. 
For Mi. and (Fy 19 the phase angles are +96 degrees and -96 degrees, 
respectively; therefore the maximum and minimum values of these components 
occur when the blade is approximately horizontal. This would be expected 


from the geometry. The phase angle for te is -159 degrees; therefore 


31a oswell, R.J., "A Method of Calculating the Spindle Torque of a 
Controllable-Pitch Propeller at Design Conditions," David Taylor Model 
Basin Report 1529 (August 1961). 
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- maximum value occurs when the blade is nearly vertical downward (6 
o'clock position). This is again as would be expected from geometry. The 
phase angles would not be expected to be precisely +90 degrees or 180 
degrees since the propeller has 22 degrees of skew. The amplitudes of 
(Fi. and uma ee each should be equal to the combined weight of the 
blade and that portion of the appropriate flexure at radii greater than 

- the radius of the appropriate strain gage. The values of ie and 
F142 were confirmed by weighing the blade and appropriate flexures. 

The values of Mis and ig are essentially zero because the blade is 
skewed and raked so that mass of the blade is balanced about the spindle 
axis in both the plane containing the spindle axis and the propeller axis, 
and the plane normal to the propeller axis which contains the spindle 

axis (see Figure 2 and Table 1). The value of E46 is nearly zero since 
F. is always in a nearly horizontal direction. For all components the 
experimentally determined amplitude and phase of the first harmonic were 
essentially independent of rotational speed n. The experimental air-spin 
results were faired so that values of the amplitude and phase of the first 
harmonic of the loading components used for separating hydrodynamic loads 
from total loads were constant, independent of n. 

Approximate scaling parameters for centrifugal loads are (F/pgn’D*) 
and (M/ppn’D”), whereas appropriate scaling parameters for gravitational 
loads are (F/pp@D°) and (M/p ped"). The model experiments presented in this 
report were conducted at full scale values of Froude number rs = (V/vgL) 


and advance coefficient J = (V/nD). Constant Froude number implies that 


V~(gL) ?~(gp)# 
ve~ gD 
Constant advance coefficient implies that 
V~nD 
v~0n 
Therefore, 
g~n-D 
ppgD°~ ppn’p* 
2n0 


P psd ~Ppn D 


Thus, for the results presented in this report gravitational loads and 


centrifugal loads scale the same. Furthermore, if proper allowance is 
made for the difference in density between the model propeller and the 
full scale propeller,* the gravitational and centrifugal loads scale the 
same as the hydrodynamic loads. 

Therefore, in addition to the components of loading arising from 
hydrodynamic effects alone, for many experimental conditions the compo- 
nents of loading arising from the sum of hydrodynamic, centrifugal, and 


gravitational effects are presented. The components of loading arising 


from the sum of the hydrodynamic, centrifugal and gravitational effects 

are designated components of total loading. The centrifugal and gravi- 
; tational loads presented are equivalent values for a nickel-aluminum- 

bronze propeller blade. These combined, or total, loading results are 


discussed in later sections. 


The time-average centrifugal spindle torque results, M, are com- 


c 
pared in Figure 9 with calculated values using the method of Roapedils 


Previous measurements of spindle torque by Boswell et ain and by Hawdon 
et al” have correlated well with values calculated by this procedure. 


The calculated value of M, is 33 percent lower than the experimental 
c 
value at design pitch (see Figure 9); however, this is within experimental 


accuracy. The largest measured value of uM, is 0.5 inch-pounds (0.07 N-m) 
(ol 
as shown in Table 6 whereas the accuracy of this measurement is plus or 


minus 0.2 inch-pounds (0.02 N-m) as discussed in the section on accuracy. 


The large experimental inaccuracy as a percent of the measured M, value 
c 
results from a combination of (1) the small value of the measured M_ F 
c 
and (2) the inferior characteristics of flexure number 3, which measures 


*The model propeller used in these experiments was made of aluminum, den- 
sity pp=5.44 lbf-s2/ft4 (2.80 g/cm3). The full scale propeller on the 
DD-963 is made of nickel-aluminum-bronze, density Pp=14.48 lbf-s2/f£t4 
(7.46 g/cm3). 

32a Oswell, R.J. et al, "Experimental Spindle Torque and Open-Water Per- 
formance of Two Skewed Controllable-Pitch Propellers," David Taylor Naval 
Ship Research and Development Center Report 4753 (December 1975). 
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Fo and M> relative to the other two flexures as discussed in the section 


on calibration. 


INFLUENCE OF DYNAMOMETER BOAT 

The results of the wake surveys with and without the downstream body 
(dynamomemter boat) are presented in Figure 10, and in Appendix A. These 
data indicate that the downstream body had only a small effect on the 
circumferential and radial variation in the flow and only a small effect 
on the harmonic content of the flow. However, they also indicate that the 
downstream body reduced the volume mean velocity through the propeller 
disk by approximately 12 percent; i.e., without the downstream body the 
volume mean wake (1-w),)=1.06 and with the downstream body (1-w),)=0.93. 
These results are, of course, without the propeller in place. 

The change in effective velocity through the propeller due to the 

downstream body was deducted from thrust and torque identities between the 
mean thrust and torque measured during the blade loading experiments at 
the self propulsion point (Condition 1 in Table 3), and mean thrust and 
torque measured during a previous self propulsion model experiment.* 
These results, which include the effect of the propeller, indicate that 
the downstream body reduced the effective velocity through the propeller 
disk by approximately 5 percent; i.e., without the body, (1-w,)=1.02 and 
(1-w))=1.00, whereas, with the body, (1-w,,)=0.97 and (iow, )=0.95. 

The difference between the effect of the downstreain body on volume 
mean wake and effective wake is probably due to a combination of the 
following: 


1. The effect of the propeller action; (1-w,) and (l-w.) include 


Q 
the effect of the propeller but (1-wyy) does not. 

2. Experimental inaccuracies; both methods for calculating the 
change in velocity are based on a small difference of two much larger 


nearly equal experimental results. 


*DTNSRDC experiments 21 and 22 on Model 5265-1B, in which the mean 
thrust and torque was measured using transmission dynamometers mounted 
inside the model hull. 
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It is judged that the dominant cause of the discrepancy is the effect of 


the propeller. 

Based on these results it was concluded that the downstream body 
reduced the mean velocity into the propeller by 5 percent at the self- 
propulsion condition. This is somewhat smaller than the 10 to 14 percent 
reduction in effective wake that was obtained in Reference 2 in which the 
same dynamometer boat was used behind a single screw model hull. It was 
assumed that the 5 percent reduction in the present expeciment occurred 
at all conditions at which experiments were conducted. Therefore, after 
the effects of centrifugal force were subtracted from the measured loading 
components as discussed previously, the time-average value per revolution 
of each hydrodynamic loading component was corrected for the effect of the 
downstream body as follows: From the measured hydrodynamic blade thrust 


(FY ) and hydrodynamic blade torque ee effective advance coefficients 
H 


based on thrust identity (Jp) and torque identity (Jy) were deduced from 

the open-water data (Figure 11). These values were multiplied by (1/0.95) 
to obtain corrected values of Jp and Jog: i.e., without the downsteam body. 
The corrected values of F_ and M_ were then obtained from the open-water 


data at the corrected advance coefficients Jn and Ja respectively. It 


was assumed that the downstream body did not affect the radial center of 


thrust F_ and tangential force FS . Therefore, 
H 


M_ corrected = (F_ corrected/F measured) eo measured) 
H H 


F corrected = (M_ corrected/M measured) (F_ measured) 
H *y y Vy 


The spindle torque (M, ) was corrected by the same procedure as used 
H 


for Lem and a using unpublished hydrodynamic spindle torque data for 


the DD-963 propeller. No corrections were made to a for the effect of 
H 
the downstream body; however, i is very small for all experimental con- 
H 
ditions, as discussed later. 


No correction for the effect of the downstream dynamometer boat was 


made to the measured circumferential variation of the loading components. 
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Calculations made by the methods of Tsakonas et a}? and McCarthy”° indi- 


cated that the influence of the downstream body alters the peak-to-peak 
circumferential variation of the loads by no more than 2 percent. How- 
ever, these methods did not agree well with the experimental results, as 


discussed in the section on correlation with full-scale data and theory. 


STEADY-AHEAD OPERATION 

For operation near the self-propulsion point (Condition 1 in Table 
3), Figure 12 presents the variation of the various components of total 
blade loading with blade angular position and Figure 13 presents the 
amplitude of the first 25 harmonics of the various components of total 
blade loading. 

Based on the dynamic calibration, as discussed in the section on 
calibration, it was judged that for all loading components the data are 
valid for the first 10 harmonics. In addition, the wake data shows no 
significant amplitudes for harmonics greater than the tenth; see Appendix 
A. Therefore, all data and analysis except Figures 12 and 13 are based 
on reconstructed signals using the first 10 harmonics. The symbols shown 
in Figure 12 indicate unfiltered values determined from the experiment; 
each represents the average value at the indicated blade angular position 
for over 200 propeller revolutions. The variation in measured values at 
a given angular position is discussed in the section on accuracy; see 
Figure 8. The lines on Figure 12 are the signals reconstructed from the 
first 10 harmonics. Figure 12 indicates that the variation of the 
signals with blade angular position are adequately represented by the 
number of harmonics retained. 

Figure 13 shows that there are no significant resonances for any of 
the loading components below the 23rd harmonic, which corresponds to 
(23)x(14.08)=324 hertz. This is higher than the lowest frequency reson- 
ance obtained in the results presented in Reference 2; i.e., 247 hertz. 
As discussed in the section on calibration, the higher frequency of the 
fundamental significant resonance obtained in the present experiment was 
anticipated because smaller and lighter blades were used in this experi- 


ment than were used in Reference 2. 
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The variation of all measured loading components with blade angular 
position for the self propulsion condition (Condition 1 in Table 3) is 
shown in Figures 14 and 15 for the hydrodynamic loads, and is shown in 
Figure 16 for the total (hydrodynamic, centrifugal and gravitational) 
loads. The amplitudes and phases of the harmonics of these loading compo- 
nents are presented in Figure 17 for the hydrodynamic loads and in Figure 
18 for the total loads. Appendix B presents tabulated values of all the 
data in Figures 14 through 18, and Table 7 presents a summary showing the 
maximum value, minimum value and first harmonic of each loading component. 
The values for each loading component are presented as decimal fractions 
of the time-average value of the corresponding loading component. These 
time-average values for both hydrodynamic loads and total loads are pre~- 
sented in Table 8. 

These data show that for hydrodynamic loading the variation of all 
loading components was predominantly a once-per-revolution variation. The 


extreme values for all loading components, except F, and My» occurred near 


the angular position of the spingle axis, 6=124 oe S degrees; i.e., 34 
degrees beyond the horizontal. At these positions the blade tip is 
approximately 12 degrees beyond the horizontal. This suggests that the 
tangential component of the wake is the primary driving force; see Figure 
10. The extreme values of Fy and M, occur at up to 120 degrees after the 
extreme values of the other components. The reason for this variation in 
location of extreme values is not clear; however, it may be partially due 


to experimental inaccuracy with the F -M, flexure as discussed in the 


section on calibration. : 
For total loading, the variation of all components with blade 
angular position follows basically the same pattern as for hydrodynamic 
loading. This occurs because the unsteady loading due to gravity, which 
is a pure first harmonic of blade angular position, is much smaller than 
that due to hydrodynamic force. Further for all components with a mea- 
surable gravitational load, except Fo the maximum value occurs near the 
angular position at which the spindle axis is horizontal; i.e., the 


gravitational load is nearly in phase or 180 degrees out of phase with 


the hydrodynamic load. 
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For hydrodynamic loading, F_ and x, were the largest measured force 
H 
and moment components; see Table 8. For these components the maximum 


values were approximately 1.43 times the time-average values, and the 
maximum value minus the minimum value (double amplitude) was approximately 


0.91 times the time-average values; see Table 7. For F. and|M the 
H 
maximum values and range of values with blade angular position were 


slightly smaller fractions of the respective time-average values. For 


x. » the maximum value and range of values with blade angular position 
H 
were much greater fractions of its time-average. This large fractional 


variation in F_ occurs because|F_ |was very small. For|M_ j|the maximum 
74 74 74 
value and range of values were 1.30 and 0.67, respectively, times the time- 


average value. The radial point of application of FY varies from 0.68R 
H 
to 0.73R, and the radial point of application of — varies from 0.67R to 
H 
0.79R. 


The maximum values of and K, were approximately 1.41 times the 
time-average values, and the range of values were approximately 0.88 times 
the time-average values for combined hydrodynamic, centrifugal and gravi- 
tational loading components, or total loading compouents. These results 
are nearly the same as the hydrodynamic results since the centrifugal and 
gravitational loads are small for these components; see Tables 6 and 8. 


The centrifugal and gravitational loads are a significant portion of 


the total loads for other loading components; for total loads Fr, (Enid 
MAX 
whereas for hydrodynamic loads F /F_=1.40. This large difference 
Yu,MAX YH 
results from the combination of centrifugal force adding to the time- 


average hydrodynamic force and gravitational force substracting from the 
circumferential variation of hydrodynamic force. Similarly, for total 


loads M /M,. =1.26 whereas for hydrodynamic loads/M /M. |=1.34. The 
*MAX MAX *Y 
centrifugal loading F, is two orders of magnitude larger than the time- 


average hydrodynamic loading, and for M,, the centrifugal loading is almost 


as large as the time-average hydrodynamic loading; see Table 8. 


The results presented here follow trends similar to those in Refer- 
ence 2 for the FF-1088 which is a single screw transom stern configura- 
tion. The component Ms which is the largest moment component for both 


cases, yields M, M, =1.40 for the present configuration (DD-963 Class) 
MAX 
and 3 /M,=1. 38 for Reference 2 (FF-1088). The maximum and minimum 
MAX 
values occur at approximately the same angular position of the blade mid- 


chord at the 70 percent radius for the two configurations. 


HULL PITCH 


Figure 19 presents the variation of the peak values and time-average 
values per revolution of the various components of blade total (hydro- 
dynamic, centrifugal and gravitation) loading* with hull pitch angle y for 
both quasi-steady simulation (time rate of change of hull pitch angle 
~=0) and unsteady simulation (}#0). These data show that for the quasi- 
steady simulation the time-average value per revolution of each loading 
component remains within 6 percent of its value corresponding to self- 
propulsion in calm water. The time-average value per revolution for the 
unsteady simulation, deviates by up to 12 percent from its value corre- 
sponding to self propulsion in calm water. 

Data at each specified value of hull pitch angle y for the quasi- 
steady runs were recorded and averaged for a minimum of 200 propeller 
revolutions whereas data for the dynamic pitching runs at each specified 
y represented an average of from 10 to 35 propeller revolutions. As dis- 
cussed earlier, the selection of a propeller revolution at a specified yp 
during the dynamic pitch runs necessitated a tolerance of only 0.05 degree 
to »y. Therefore, the differences between the results for the quasi-steady 
and unsteady simulations, including the time-average values per revolution, 
were significantly larger than any errors which may have arisen from inac- 
curacies in setting the experimental conditions. 

For quasi-steady simulation, the absolute value of the time-average 


value per revolution of all loading components, except spindle torque M> 


* 

No results are ahown for F_ since the F_ loading arises primarily from 
centrifugal effects, as discussed previously, which are independent of 
hull pitch. 
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were larger for the stern-up condition than for the stern-down condition; 


the largest value occurs at (yp )=1.85 degrees or time=0.31 seconds in 


- 
the reference of Figure 19. 2D og nal that the effective speed of 
advance of the propeller increases slightly for the stern-down condition 
and decreases slightly for the stern-up condition. This appears reason- 
able since for stern-up the propeller tends to be further into the bound-~ 
ary layer of the hull. However, the time-average value per revolution did 
not monotonically increase with increasing » for all components. 

For dynamic simulation the largest absolute value of the time- 
average value per revolution of ali loading components, except spindle 
torque My, occurs at approximately 0.15 second after the condition 
(WY =0 ~>0, which is the reference for time t=0 in Figure 19. This 
indicates that the maximum time-average value during dynamic simulation 
occurs at a value of hull pitch angle » which occurs 0.16 second or 0.1 
cycle, before the p at which the maximum time-average value occurs during 
quasi-steady simulation. 

There was a significant difference between the peak values for the 
quasi-steady simulation and the unsteady simulation. For the quasi-steady 
simulation, the variation of the peak values with hull pitch angle yp 
followed approximately the same trends as the variation of time-average 
values per revolution. These quasi-steady results indicated that for 
Vow up to 1.85 degrees, the maximum increase in the peak value of any 
loading component above the corresponding value for v=VoW was 5 percent. 
For the dynamic simulation, however, the maximum value of the peak loads 
increased as much as 23 percent above the corresponding value for steady 
ahead at a fixed hull pitch v= Vow" 

The dynamic simulation exhibited a dramatically different trend of 
peak load with y than was indicated by the quasi-steady simulation. For 
the dynamic simulation, the largest value of the peak loading, for all 


components except spindle torque M,, occurred at approximately time t=0.8 


Z 
second in the hull pitch cycle shown in Figure 19. This corresponds to 


w=1.5 degrees stern down during the portion of the cycle in which the 


stern is moving down; i.e., (-Y oy =-1-5 degrees, ¥<0. For dynamic 


simulation, the smallest value of peak loading occurred near v= Voy as the 
hull passed from the stern-down to the stern-up portion of the cycle; i.e., 
(W-Vo)=0, ¥>0. 

This difference in the unsteady loading between the quasi-steady and 
unsteady simulations may be due to an additional relative velocity compo- 
nent arising from the motion of the hull during dynamic pitching. As the 
hull passes through v= Vow? the vertical velocity of the hull (and propel- 
ler) is a maximum. As the hull goes from stern up to stern down through 
v=Voy> the upward velocity component relative to the propeller in the 
plane of the propeller tends to increase above the values at fixed hull 
pitch at Y= Vous This tends to increase the amplitude of the first harmo- 
nic of the tangential velocity, and thereby increase the unsteady loading 
(and increase the peak loading). The maximum vertical velocity of the 
propeller for sinusoidal pitching with (ax cw 2+ 85 degrees and fre- 
quency=0.8 hertz is approximately 1.47 feet per second (0.448 m/s). This 
is equivalent to an additional tangential velocity ratio (v/v) of 0.133. 
For p fixed at V=Voy? ((V,),/V)=0.130 (see Appendix A). Therefore, 


CV) /Y max, #0 = 0.130 + 0.133 - 2.02% 


((V,) 0.130 


U”) 0,0 Voy 
This maximum occurs at a model simulated time of approximately 0.2 second 
before the maximum measured loads. The measured increase in unsteady loads 
arising from dynamic pitching was somewhat smaller than this calculated 


increase in tangential velocity, for example: 


F . = m 
MAX , #0 v= Vow = 0.6 


_ 


2 
F ue 0.44 
"MAX,V=0 “Y= Woy, 


= 1.41 


*A numerical error was found in a similar calculation presented in Ref- 
erence 2. With the numerical error corrected the results in Reference 2 
are substantially the same as those presented here. 
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On the basis of two-dimensional quasi-steady theory, the increase unsteady 
loading should be approximately proportional to the increase in tangential 
velocity.* 

The unsteady loading is important from consideration of fatigue of 
the propeller blades and hub mechanism. Since a ship may operate for an 
extended period in a seaway, the effect of the ship motions, such as 
dynamic hull pitching, on unsteady blade loads is significant. The differ- 
ence between the peak load and the time-average load per revolution is a 
measure of the unsteady loading. With this difference as a measure of the 
unsteady loading, the quasi-steady simulation indicates that for hull 
pitch angles v-Vow up to 1.85 degrees, the unsteady loading for Ms which 
is the largest moment component, increased by 8 percent above its corres- 
ponding value for vVoye By contrast, the dynamic simulation showed the 
unsteady loading for the x, component increased by 50 percent above its ; 
corresponding value for Vow without hull pitching. This indicates that 
the quasi-steady simulation is completely inadequate for estimating the 
effect of the seaway on unsteady loading. This also shows that the effect 
of the ship motions can dramatically increase the unsteady loading on the 
blades. Therefore, the effect of the ship motions due to operation in a 
seaway should be considered in any analysis of blade loading and in any 
fatigue analysis of the propeller blades or hub mechanism. 

The results presented here for hull pitching generally agree with 
the same type of results presented in Reference 2 for a model of the FF- 
1088, which is a single screw transom stern configuration. For M , which 
is the largest measured moment component in both cases, the comparative 
results, presented as a fraction of the time-average value without hull 


pitching, are as follows: 


*This simple analysis provides an upper bound to the dynamic pitching 
load, since the hull boundary above the propeller would tend to reduce 
the dynamic pitching induced upward tangential velocity relative to the 
propeller. 
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DD-963 FF-1088 

(Present (Reference 

Report) 2) 
Peak value, #0 1.60 1.57 
Peak value, j=0 1.43 1.40 
Peak value without pitching  =0, vVow 1.40 1.36 
Maximum time-average value, W#0 1.10 1.03 
Maximum time-average value, W=0 1.05 1.04 


The variation of the loading components with simulated time during the 
pitch cycle are also somewhat similar for these two configurations. The 
comparative results for Ms presented as time in seconds following the 


point V-Vqy7> W>0 are as follows: 


DD-963 FF-1088 

(Present (Reference 

Report) 2) 
Peak value, W#0 0.77 0.62 
Peak value, $=0 0.31 0.31 
Maximum time-average values, #0 0.20 0.72 
Maximum time-average value, j=0 0.31 0.31 


ACCELERATION 

The variation of all measured loading components with blade angular 
position for the quasi-steady simulated acceleration condition V=P=n=0 
(Conditions 7 through 11 in Table 3) is shown in Figures 20 and 21 for the 
hydrodynamic loads, and is shown in Figure 22 for the total (hydrodynamic, 
centrifugal, and gravitational) loads. The amplitudes and phases of the 
harmonics of these loading components are presented in Figure 23 for 
the hydrodynamic loads and in Figure 24 for the total loads. Appendix B 
presents tabulated values of the data in Figures 20 through 24. The 
values for each loading component are presented as decimal fractions of 
the time-average value of the corresponding loading component at the self 
propulsion condition (Condition 1 in Table 3). These average values for 


both hydrodynamic loads and total loads are presented in Table 8. 
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Figure 25 presents the Taylor wake fraction based on thrust 1-w,, and 
the Taylor wake fraction based on torque 1-Wo as derived from the measured 


values of ye and M. and the open-water characteristics of the propeller 


(Figure 11). These data indicate that (1-w,) varies by only approximately 
3 percent during the simulated acceleration. The value of (-w9) varies 
by only 1 percent for simulated time t>40 seconds; however, the value of 
(1-w9) varies substantially during the initial portion of the simulated 
acceleration (t<40 seconds). 

Figures 20 and 22 show that for all measured hydrodynamic and total 


loading components, except Fy which is small, the peak values, including 
H 
variation with blade angular position occurred at the self propulsion con- 


dition. That is, for the acceleration condition simulated (see Figure 7 
and Table 3), the propeller is not exposed to higher peak loads than those 
to which it is exposed during full-power steady-ahead operation. 

In contrast to the peak loads, for most components the largest time- 
average loads per revolution occurred at the first experimental condition 
(V=2.65 knots, n=10.21 revolutions per second, J=0.64) during the simulated 
acceleration maneuver. The largest measured hydrodynamic force and moment 


components, F_ andM_, yield (F /¥ )=1.21 and (M /M )= 
Hy YH ,MAX “H,SP Yu,MAX ‘H,SP 
1.29, whereas for total loading (F. /F_-)=1.16 and (M /M  )=1.21. 
*sp Ymax sp 
The conditions (F /¥ )> FF /F +) and (M /M )> 
,MAX ™“H,SP *max *sp Yu,MAX YH,SP 
(M /M_ ) occur because the centrifugal and hydrodynamic components are 
Ymax ‘sp 
additive for FY and a, (i.e., they have the same signs) and the hydro- 


dynamic loads increase with decreasing rotational speed n, whereas the 
centrifugal loads decrease with decreasing n. 
Higher time-average and peak loads than those shown in Figures 20 
and 22 could, of course, be developed during acceleration maneuvers, j 
depending on values of v, n, and P. 
Figure 21 shows the variation in the radial center of longitudinal 


and radial center of tangential force, Tr These results 


*y Yu 


show that the time-average radial centers of these force components vary 


force, Tp 
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monotonically with advance coefficient over the range evaluated. As the 


advance coefficient based on thrust effective wake, =J (1-w,) increases 


J 

i ea. 

from 0.63 (at V=2.65 knots) to its design value of 1.14 (at V=6.52 knots), 

Te decreases from 0.76R to 0.71R whereas Tp increases from 0.66R to 
*y Yu 

0.74R. ‘This variation in ¥F is the reason that (FY. /¥F )< 


a H,MAX “H,SP 


(M /M ) as discussed in the preceeding paragraph. 
Yu,MAX H,SP 


For all loading components, the variation with blade angular posi- 
tion tended to be dominated by the first harmonic for all conditions 
throughout the simulated acceleration maneuver. For all conditions at 
which there was significant variation in loading with blade angular posi- 
tion, the maximum and minimum values for all components except eS and M, 
occurred for the blade spindle axis near 6=135 or 315 degrees (blade tip 
near 6=115 or 195 degrees). This suggests that the variation in loading 
with blade angular position is produced primarily by the circumferential 
variation of the tangential velocity in the propeller plane (see Figure 
10). The angular variation of each loading component retained basically 
the same shape independent of speed and advance coefficient. 

There was a dramatic reduction in the circumferential variation of 
all measured loading components with decreasing speed V and decreasing 
rotational speed n. Previous data have shown that for a given propeller 
in a given flow field, the circumferential variation in the hydrodynamic 
loading varies approximately as the product of ship speed V and rotational 
speed n; see Wereldsma. >” Figure 26 presents results in a form which 
allows evaluation of how closely the measured unsteady loading varies with 
nV. The ordinate is the first harmonic of the components of hydrodynamic 


blade loading except Fe » which is very small, and the abscissa is nV. 
H 
The data shown in Figure 26 indicate that the first harmonic of each of 


the presented hydrodynamic loading components is approximately proportion- 
al to nV. 


s3yereldsma, R., "Tendencies of Marine Propeller Shaft Excitation," 


International Shipbuilding Progress, Vol. 19, No. 218, pp 328-332 (Octo- 
ber 1972). 


Figure 27 presents the variation of the time-average values per 


revolution and peak valves of the various components of total blade 
loading for both quasi-steady simulated acceleration (V=n=P=0) and unsteady 
simulated acceleration (V>0, Pn 

There was only a small variation in the measured loading components 
between the quasi-steady simulated acceleration and the unsteady simulated 


acceleration. For all components except M_, the largest variation between 


’ 
the results from the two types of sieaiaed oa expressed as a decimal fraction 
of the corresponding time-average value at the self-propulsion point was 
0.05 for the peak values and 0.02 for the time-average value per revolu- 
tion. The corresponding variations for M, were no greater than 0.06 for 
the peak values and 0.05 for the time-average value per revolution. 

The variation in the results between the two types of simulation 
appeared to be essentially random. This suggests that these deviations 
are some measure of the experimental accuracy and do not represent any 
systematic trends arising from the difference in V between the two types 
of simulation. 


Data for the quasi-steady simulation were recorded and averaged for 


a minimum of 200 propeller revolutions, whereas data presented for the 
unsteady runs represent an average of only five revolutions. Further, 
the steady experimental conditions which were set during the quasi-steady 
simulation allow the values of V and n to be controlled more precisely 
than during the unsteady runs; however, the average of the five values of 
V and n during the unsteady runs for which data are presented was general- 
ly within one percent of the target values. 

The results presented in this section for a simulated acceleration 
maneuver follow trends similar to those in Reference 2 for a simulated 


crash-forward maneuver on a model of the FF-1088. Both sets of data show 


the following: 


1. The values of (1-w,.) and (i-w9) do not vary substantially except 
during the initial stages of the acceleration or crash-forward maneuver. 
2. The variation of all loading components with blade angular posi- 


tion was dominated by the first harmonic throughout the simulated maneuver. 


3. The amplitude of the first harmonic of all loading components 
varied essentially as nV. 

4. The acceleration of the hull did not have a significant effect 
on the loads; i.e., the loads for quasi-steady acceleration V=0 and un- 
steady acceleration v>0 were not significantly different. 

The ratios of either the peak or time-average loads during the accelera- 
tion (or crash-forward) maneuver to the time-average loads at the self- 
propulsion point do not agree closely for the results in the present 
report (DD-963 Class) and Reference 2 (FF-1088). This difference is to 
be expected since these ratios are very sensitive to the value of v, n, 
and P for the simulated maneuvers, which are quite different for these 
two cases. The largest moment component for both cases M , gives 


/M = 1.40 from the present report and 1.51 from Reference 2. The 
Ypeak Ysp 
“results from the present report show M /M. = 1.21 whereas the results 
is . MAX ‘SP 
from Reference 2 yield M /M. =ipk 39. 
Ymax SP 
CORRELATION WITH FULL-SCALE DATA AND THEORY 

For operation near the self-propulsion point (Condition 1 in Table 
3), correlation was made between the model experimental loads obtained in 
the present investigation, bending moments deduced from strains measured 
on the corresponding full-scale propeller, and analytical calculations. 

These comparisons were made for My 3 and Mo A which were calculated 
from the thvee measured forces and three measured moments. As discussed 
in the section on experimental results, My 3 is defined as the bending 
moment applied on the spindle axis about the axis intersecting the spindle 
axis at r=0.3R and parallel to the expanded pitch line at r=0.3R. The 
vector as defined by the conventional right-hand rule for a right- 


M 
0.3 
hand propeller (left-hand rule for a left-hand propeller) intersects the 


plane normal to the propeller axis at the angle bo — tan" (P, 3/ (0. 3mD)) 


and is directed so that a positive My 3 puts the face (pressure side) of 


the blade in tension. The component M is defined in a manner analogous 


0.4 
to My 3? except it, is referred to r=0.4R. 


The full-scale data used for correlation are preliminary values of 
strains ¢ measured* at several chordwise stations at r=0.35R and r=0.45R 
on both sides of the blade on the DD-963 CP propeller. Both time-average 
strains per revolution and variation of strain with blade angular posi- 
tion were recorded during full-power stready-ahead operation. By inter- 
polation, values of radial strain at the spindle axis at r=0.4R were ob- 
tained. Assuming that the variation in radial strain is proportional to 
the variation in total (hydrodynamic, centrifugal, and gravitational) 


bending moment; i.e., that (e /€ )= (M /M. ,), these full- 
10. 4 yay r0.4 0.4 sax 0.4 
scale data indicate that (M /M. ,) = 1.48. From the model data 
0.4 vax 0.4 


(My. ae 3 


of the variation of the full scale strain with blade angular position indi- 


) = 1.43 and (M /M, ,) = 1.56. A cursory evaluation 
0.4 ax 0.4 


cates that it follows trends similar to the bending moment determined from 
the model experiments. The correlation with full scale data presented 
here is preliminary; a more thorough correlation with the full scale data 
will be undertaken when analysis of the full scale data is complete. 
Theoretical calculations of hydrodynamic loads were made by using 


the method of Tsakonas et “aoe 


which is based on unsteady lifting- 
surface theory and the method of MeCaitthy , 7° a quasi-steady technique 
which utilizes the open-water characteristics of the propeller. Although 
the method of McCarthy is a simple quasi-steady technique, it was judged 
that this method should be suitable to the cases under consideration in 
this report because the dominant unsteady loading occurs at a low reduced 


frequency and the dominant first harmonic of the wake is in phase radially 


*The full-scale measurements were conducted by personnel in DTNSRDC 
Code 1962 under the direction of G.P. Antonides. The results presented 
here are preliminary, and thorough analysis of the data is continuing. 
The details of this full-scale trial will be reported in a future DTNSRDC 
report. 


*"eeakonas, S. et al, "Documentation of a Computer Program for the 
Pressure Distribution, Forces and Moments on Ship Propellers in Hull 
Wakes," (In Four Volumes), Stevens Institute of Technology, Davidson 
Laboratory Report SIT-DL-76-1863 (January 1976). Revised April 1977. 
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(see Figure 10 and Appendix A). These calculations were made for Condi- 


tion 1 in Table 3 using the wake measured in the plane of the propeller 
both with and without the downstream dynamometer boat in place (Figure 10 
and Appendix A), and with the mean velocity through the propeller deter- 
mined from thrust identity used as the reference velocity. 

The use of the speed of advance based on thrust effective wake, 
V,=V(1-w,), as the reference speed in these calculations is consistent 
with the use of this velocity to correct the time-average loads for the 
effect of the dynamometer boat as discussed in the section on experimental 
results. Tsakonas et ao" recommend using the ship speed as the refer- 
ence velocity, which is equivalent to using (1-w,)=1.05 however, this 
recommendation was not followed here because the flow does not pass through 
the propeller at the ship speed. To evaluate the sensitivity of the pro- 


cedure of Tsaknoas et al,7°">" 


to the reference speed, calculations were 
performed for the first harmonic using the thrust effective wake (1-w,,) 
and the volume mean wake determined from the wake surveys (1-wy)) + These 
calculations showed the following: 

Wake without dynamometer boat 


M (using (l-w,,) = 1.02) 
va ss = 0.99 
(using (1-wy) = 1.06) 


M 
0.3) 


*yo.3, Or? = *yo. 3, “vw = -0.3 degrees 


Wake with dynamometer boat 


“0.3, (using (1-w,) = 0.97) 


Mo. 9 (using (1-4) = 0.93) 


= 1.02 


Sho. 3 Mr ' *yo.3, “ww = 0.5 degrees 


Therefore, the calculated unsteady loads using the method of Tsakonas et 


ac are not sensitive to the reference speed over the range of con- 


cern in the present case. 
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For the method of Tsakonas et al, calculations were conducted 
for the first ten harmonics of the wake,” The "normal" components of wake 
harmonics, as required by this method, were defined as the wake harmonics 
normal to the chord line of the blade section at the local radius rather 
than normal to the advance angle at the local radius as recommended by 


Tsakonas et al.??9 


With the wake harmonics resolved normal to the blade 
chord, this method apparently considers both the unsteady flow normal to 
the resultant inflow and an approximation to the unsteady flow parallel to 
the resultant inflow. 

The quasi-steady calculations are based on the circumferential vari- 
ation of the wakes measured at the 0.7 radial station. These calculations 
were made at 10-degree increments of blade angular position 8. It is 
assumed that the radial centers of the unsteady thrust and tangential 
force are at r/R=0.70 for all blade angular positions. 


Figure 28 presents values of M and M with and without the 
0.3, 0.4, 


downstream dynamometer boat, calculated with the methods of Tsakonas et 
ai??? and McCarthy. 2° Based on these calculated results it appears that 
the dynamometer boat does not have a significant influence on the cir- 
cumferential variation of the blade loads. 

Figures 29 and 30 present the variation with blade angular position 


and the first ten harmonics of My 3 and My 4 from the model experiments 
34 45 


and analytical calculations. All data are nondimensionalized on the same 
quantity, i.e., the time-average bending moment determined from the model 
experiments and corrected for the downstream body as discussed in the sec- 
tion on experimental results. This comparison indicated that the experi- 
mental unsteady hydrodynamic bending moments were substantially higher than 


the calculated results. A typical comparison is as follows: 


*These calculations were made by using the computer program developed by 

the Davidson Laboratory including refinements made through April 1977. 
None of the refinements made since December 1975 influence the calculated 
unsteady loads presented in this report. Therefore, this calculation pro- 
cedure is the same as that used in Reference 2 for calculating the unsteady 
bending moments on the propeller on the FF-1088. 

> nsakonas, S. et al, “Correlation and Application of an Unsteady Flow 
Theory for Propeller Forces," Transactions of the Society of Naval Archi- 
tects and Marine Engineers, Vol. 75, pp 158-193 (1967). 
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Prediction Method 


M -¥M M -M 
0.3, MAX 0.3, 0.4, MAX 0.4, 
M3 (experiment) Moy (experiment ) 

H H 
Model Experiment 0. 36 0.36 
Quasi-Steady 
Procedure26 0.25 0.24 
Unsteady Procedure” 0.20 0.19 


For this typical comparison the experimental result is approximately 47 
percent higher than the quasi-steady prediction and approximately 85 per- 
cent higher than the unsteady prediction. 

The circumferential variations in the model experimental results of 


other components of blade loading F_. , F_,M_, and M_ were larger than 
a: A uy 


the values calculated by the two indicated procedures by approximately 


the same ratio as shown for se and MO. 


These comparisons are not 
shown. 

Previous investigators have compared experimental unsteady forces 
and moments on a single blade of various propellers in inclined flow with 
forces and moments calculated by a quasi-steady procedure similar to that 
described by McCarthy. 7° These. experimental loads were obtained by direct 
measurement of unsteady forces and moments on a single blade (References 
2, 3, 17, 18, and 19) or were deduced from measured steady transverse 
forces and moments along axes fixed relative to the flow, i.e., not rotat- 
ing with the propeller (Reference 36). References 2, 3, 17, 18, 19, and 
36 all show that for noncavitating conditions, the experimental unsteady 
blade loading was from 1.5 to 2.0 times as large as the values calculated 
by the quasi-steady method. This agrees with the results of the present 
investigation; see Figures 29 and 30. 


sovesche, F., "The Study of Ships' Propellers in Oblique Flow," 
Defence Research Information Centre Translation No. 4306, Copyright Con- 
troller: Her Majesties Stationary Office, London, England, October 1975; 
English Translation of "Untersuchung von Schiffsschrauben in schrager 
Anstromung," Schiffbauforschung, Vol. 3, No. 3/4, pp 97-102 (1964). 
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Preliminary results from blade loading experiments”” conducted in 
idealized flows indicate that: 

1. The unsteady blade loads in either axial or tangential wakes are 
not influenced by the presence of a nearby boundary above the propeller. 

2. In inclined flow without an upstream hull, the experimental un- 
steady loading near the design advance coefficient is nearly two times as 
large as that calculated by the method of Tsakonas et ai, 2934 and approxi- 
mately 80 percent larger than that calculated by the method of McCarthy. 7° 

3. In an axial wake with a dominant first harmonic of blade angular 
position which was generated by upstream wire grid screens with zero shaft 
angle, the unsteady loading near the design advance coefficient is within 
approximately 1s ee of the values calculated by the methods of 

’ 


Tsakonas et al, and Mecaeeny.?” This is in agreement with previous 


correlations with the method of Tsakonas et al for unsteady shaft (bearing) 
forces and moments for operation in axial wakes, °?12 
These results indicate that the large discrepancy between the experi- 
mental and calculated unsteady bending moments presented in the current 
report appear to be due to the inability of the present theories to 
properly account for all important characteristics of the flow field for 
operation in inclined flow. All available procedures, including the un- 


steady theory of Tsakonas et ai??? 


and the quasi-steady procedure of 
MeCarthy,”° implicitly assumed that the propeller slipstream follows the 
propeller axis rather than the direction of the effective velocity into 
the propeller which is at an angle to the propeller axis in inclined row. 
It is speculated that this failure to consider the true direction of the 


slipstream is the major factor in the analytical underprediction of the 


*The results in Reference 12 were at substantially higher reduced fre- 
quency than the results in the present study; therefore, the quasi-steady 
procedure of McCarthy over-predicted the unsteady loads in Reference 12. 
3) Sewell, R.J. and S.D. Jessup, "Experimental Determination of Period- 
ic Propeller Blade Loads in a Towing Tank," Presented to the 18th American 
Towing Tank Conference, U.S. Naval Academy, Annapolis, Maryland (August 
1977). 
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unsteady blade loads in inclined flow. Numerical computations to check 


this hypothesis are planned. 


SUMMARY AND CONCLUSIONS 

Experiments were described in which the mean and unsteady loads, 
including hydrodynamic, centrifugal, and gravitational loads, were mea- 
sured on a single blade of a model of a CP propeller on the DD-963 Class 
Destroyer. The experiments were conducted behind a model of the DD-963 
hull under steady-ahead operation, hull pitching motions, and simulated 
acceleration maneuvers. The discussion of experimental techniques in- 
cludes a description of the dynamometer and data analysis system. The 
results are summarized as follows: 

1. For all significant loading components, except for radial force, 
the loads are predominantly of hydrodynamic origin. 

2. The circumferential variations of all measured components of 
hydrodynamic and total blade loading are primarily a first harmonic, with 


maximum and minimum values occurring near the blade angular position which 


is 25 degrees past the position at which a radial line from the propeller 


axis to the tip is horizontal. 
3. For steady-ahead operation: 

a. The maximum values and peak-to-peak circumferential varia- 
tions for measured hydrodynamic forces and bending moments were up 
to approximately 1.43 and 0.91 of the time-average values, 
respectively. 

b. The maximum values and the peak-to-peak circumferential 
variations for measured total forces and bending moments were up to 
approximately 1.41 and 0.88 of the time-average values, respectively. 

c. The model results for circumferential variation of bending 
moments about the nose-tail lines of the 0.3 and 0.4 radii agreed 
fairly well with loads deduced from strain measurements on the full- 
scale propeller, but they were larger than theoretically calculated 


values. 


4. For simulated hull pitching (maximum pitch angle of 1.85 
degrees): 

a. The maximum values of measured total forces and bending 
moments increased over the corresponding values without hull pitch 
by 5 percent for quasi-steady simulation and by 23 percent for un- 
steady simulat‘on with model pitching frequency equal to 0.8 hertz 
(full scale equivalent frequency is 0.16 hertz). 

b. The peak-to-peak circumferential variation of the measured 
total forces and bending moments increased over the corresponding 
values without hull pitch by approximately 5 percent for quasi- 
steady simulation and by approximately 50 percent for unsteady 
simulation with model pitching frequency equal to 0.8 hertz. There- 
fore, any quasi-steady simulation of ship motions is completely in- 
adequate for estimating the effect of ship motions on unsteady pro- 
peller blade loading. 

5. For the simulated acceleration maneuver: 

a. The dominant first harmonic of the measured hydrodynamic 
forces and bending moments varied in a nearly linear manner with 
the product of ship speed and propeller rotational speed. 

b. The acceleration of the hull did not have a significant ef- 
fect on the measured loads. Therefore, propeller blade loading 
during an acceleration maneuver can be adequately estimated by 
quasi-steady experiments. 

c. The maximum time-average values of measured forces and bend- 
ing moments per revolution, were in the range of 1.21 to 1.29 of the 


time-average values during full-power steady-ahead operation for 


hydrodynamic loads, and in the range 1.16 to 1.21 of the time-average 


values during full-power steady-ahead operation for total loads. 

d. The simulated acceleration condition did not expose the pro- 
peller to higher peak loads than those to which it is exposed during 
full power steady-ahead operation. However, these loads are very 
sensitive to the maneuver simulated and substantially higher peak 


loads could be developed during other acceleration maneuvers. 


e. Except for the initial portion of the simulated accelera-~ 
tion maneuver, the Taylor wake fractions were within three percent 
of their values at the self propulsion point. 

All of the results presented here on a model of the DD-963 Class 
Destroyer follow close to previously reported results of similar experi- 


ments on a model of the FF-1088. 
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APPENDIX A 


DETAILS OF WAKES* 


Tables 9 and 10 present the velocity component ratios and the 
harmonic content of the wakes in the plane of the propeller, both with 
and without the downstream dynamometer boat. The data at even radial sta- 
tions were obtained by interpolation and extrapolation of the measured 


data as described in.Reference 38. 


*Al1l data presented in Appendix A were obtained from wake surveys con- 
ducted by R.F. Roddy, DINSRDC Code 1524. Further details of these wake 
surveys will be presented in a future DTNSRDC report. 


38cheng, H.M., “Analysis of Wake Survey of Ship Models - Computer 


Program AML Problem No. 840-219F," David Taylor Model Basin Report 1804, 
March 1964. 
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APPENDIX B 


DETAILED EXPERIMENTAL RESULTS 


Table 11 presents detailed experimental results, including variation 
with blade angular position and harmonic analyses, for steady-ahead opera- 
tion at V=6.52 knots, n=14.08 rev/sec. The data in Table 11 are tabulated 
values of the data presented in Figures 12, 13, 14, 16, 17, and 18. 

Tables 12 to 15 present detailed experimental results including 
variation with blade angular position and harmonic analyses, for the quasi- 


steady acceleration conditions. The data in Tables 12 to 15 are tabulated 


values of the data presented in Figures 20, 22, 23, and 24. 
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Figure 5 - Typical Strain-Gaged Flexure 
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Figure 6 - Arrangement of Typical Flexure in Hub 
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Figure 7 - Experimental Acceleration Conditions 
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Figure 8 - Experimental Data Showing Plus and Minus Two 
Standard Deviations on Measured Values of F, 
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Figure 10 (Continued) 
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Figure 10b - Harmonic Content at r/R = 0.80 
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Figure 10 (Continued) 
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Figure 12 - Influence of Extraneous Signals on Measured Loads 
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Figure 12 (Continued) 
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Figure 12 (Continued) 
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Figure 12 (Continued) 
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Figure 12 (Continued) 
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Figure 13 - Experimental Data Showing Extraneous Higher Harmonics 
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Figure 13 (Continued) 
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Figure 13 (Continued) 
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Figure 13 (Continued) 
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Figure 13 (Continued) 
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Figure 20 - Variation of Experimental Hydrodynamic Loads with 
Angular Position for Quasi-Steady Acceleration 
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Figure 20 (Continued) 
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Figure 20 (Continued) 
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Figure 20 (Continued) 
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Figure 22 - Variation of Experimental Total Loads with Angular 
Position for Quasi-Steady Acceleration 


V= 2.65 KNOTS 
n = 10.21 REV/SEC 


V= 3.55 KNOTS 
n = 10.96 REV/SEC 


Fil noe 


V= 5.36 KNOTS 
n = 12.70 REV/SEC 


V= 6.26 KNOTS 
n = 13.78 REV/SEC 


V= 6.52 KNOTS 
n = 14.08 REV/SEC 


POSITION ANGLE @ (DEGREES) 


Figure 22a - fo 


109 


Figure 22 (Continued) 
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Figure 22 (Continued) 
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Figure 22 (Continued) 
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Figure 22 (Continued) 


V= 2.65 KNOTS 
n = 10.21 REV/SEC 


V= 3.55 KNOTS 
n = 10.96 REV/SEC 


Mo.3/Mo.35p 


Vv 
n 


= 6.26 KNOTS 
= 13.78 REV/SEC 


POSITION ANGLE @ (DEGREES) 


Figure 22g - Mo 3 


“0 45 90 135 180 225 270 315 360 


ri aie, 


Nel Seo coe Ne rere 


Pe peat aie ed eee abadiedbon apcesieertvigy ancbaten eae naa aa aT! ee 


Figure 22 (Continued) 


V= 2.65 KNOTS 
n = 10.21 REV/SEC 


INNES ae ne 
PSPS 
JAZ | LN OAL 
Pt | sme RS LL 
ate yee. LRU 
P| | | esate AA 


0 45 90 135 180 225 270 315 360 
POSITION ANGLE 6 (DEGREES) 


| 
. 
| 
| 
| 


Mo.4/Mo.45p 
— 
f=) 


Figure 22h - Mo 4 


116 


Figure 23 - Harmonic Content of Experimental Hydrodynamic 
Loads for Quasi-Steady Acceleration 
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Figure 23 (Continued) 
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Figure 2 3 (Continued) 
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Figure 23 (Continued) 
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Figure 24 - Harmonic Content of Experimental Total Loads 
for Quasi-Steady Acceleration 
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Figure 24 (Continued) 
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Figure 25 - Taylor Wake Fractions during 
Simulated Acceleration Maneuvers 


Figure 26 - Variation of First Harmonic of Experimental Hydrodynamic 
Loads with nV for Quasi-Steady Acceleration 
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Figure 26 (Continued) 
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Figure 27 - Comparison of Time-Average Values per Revolution and Peak Values of 
Various Components of Experimental Total Blade Loading for Quasi-Steady 
and Unsteady Simulated Acceleration 
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Figure 27f - My 3 


Figure 28 - Variation of Hydrodynamic Bending Moment at 30 Percent 
and 40 Percent Radii with Blade Angular Position, Theoretical 
Prediction With and Without Dynamometer Boat 
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Figure 28 (Continued) 
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Figure 28b - Quasi-Steady Theory 
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Figure 30 - Harmonic Content of Hydrodynamic Bending Moment at 30 Percent 
and 40 Percent Radii, Comparison of Model Data and Theory 
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Figure 30 (Continued) 
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TABLE 1 — CHARACTERISTICS OF PROPELLERS ON DD-963 CLASS 
DESTROYER; DINSRDC MODEL PROPELLERS 4660 AND 4661 


OA NORE TENOR Pee ey I TR Ne 


Diameter, D: 17.0 feet (5.577 m)* Speed at Full Power, V: 32.5 knots 
Rotation: Inward Turning* Hub-Diameter Ratio, D, /D: 0.30 
Number of Blades, Z: 5 Expanded Area Ratio: 0.73 
Maximum Rotational Speed (Rated), n: Blade Thickness Fraction: 0.054 
168 r/min Section Meanline: NACA a= 0.8 
Full Power (Rated), Pp : 
40,000 hp (29,800 kW) 
c/R c/D P/O 02° (deg) «= (Z,, 0°" vd fylc 
0.3 0.178 1.165 2.985 -0.0006 0.0420 ti) 
0.35 0.210 1.296 3.481 -0.0022 0.0372 0.0050 
0.45 0.271 1.480 4.810 -0.0095 0.0290 0.0209 } 
0.55 0.327 1.566 6.631 -0.0185 0.0226 0.0267 4 
0.65 0.374 1.566 8.978 -0.0288 0.0178 0.0256 
0.75 0.406 1.498 11.895 -0.0392 0.0146 0.0209 
0.85 0.409 1.381 15.410 -0.0488 0.0122 0.0151 
0.90 0.387 1.306 17.403 ~0.0529 0.0110 0.0122 
0.95 0.326 1.222 19.557 -0.0561 0.0091 0.0094 
1.0 0 1.128 21.876 -0.0583 0 0 


*For model propeller, D = 0.6848 feet (22.47 cm) 


“The blade loads were measured on the starboard propeller (left-hand rotation); 
DTNSRDC Model Propeller 4661 


**The spindle axis is the blade reference line. 
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TABLE 2 - CALIBRATION MATRIX 


~3.33 -0.099 0.058 -0002 -0.090 -0.066 
~0.008 5.00 -0.002 0.084 -0.168 -0.196 


-0.071 0.048 -334 0.070 -0.014 0.058 
Calibration Matrix = (C, ) = 
i 0.009 ©«-0.033 0.103 499 -0.071 0.297 
-0.234 -0.367 0.291 -0861 -999 -1.87 


where, 
F Fra Fx Fy, F., are indicated forces in votts 
My, My, M,,-M,,.M,, are indicated moments in volts 
Fy, E Fi, ee Fra. Fra’ Fra are applied forces in pounds 
M,, My, a My, My,° M,, are applied moments in inch-pound 
F., Fl, C,; for j=1,3,5 are in volts per pound 
M,, mM, C,, for j = 2, 4, 6 are in volts per inch-pound 


The arrays are arranged by the coupling of the force and moment peirs of the three flexures; 
i.e., Flexure 1 mesures F, and M,, Flexure 2 measures F, and M,, and Flexure 3 measures 


F, and M,. 
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TABLE 3 - MODEL EXPERIMENTAL CONDITIONS 


Condition n nhs Yow 
Number } knot tis Oia. 


Er CE 


oo sale 


Hull Pitch 
Unsteady Hull 
Ea CR SC 


Quasi-Steady 
Acceleration 


2.65°* (1.47) ** 
3.55** (1.96) ** 
5.36** (2.97) ** 
6.26** (3.47)** 
6.52** (3.61) ** 


(0.31) ** 
0.56** (0.31)** 
0.32** (0.18)** 
0.08** (0.04)** 
(0)** 


oo ooicoocd o& 


°o 


“Sinusoidal with amplitude equal to 1.85 degrees, frequency equal to 0.8 Hz. 
** Varies with time (Figure 7); value shown is at time of interest. 


TABLE 4 - FULL-SCALE CONDITIONS SIMULATED BY MODEL EXPERIMENTS 


Condition . Yow 
—— piste Ls Seth 


[Self Propulsion] 1 Pe ee 


Quasi-Steady 2 80 (18.0) 
Hull Pitch 32.5 (18.0) 
32.5 (18.0) 
32.5 (18.0) 


0.66** (0.31)** 
0.56** (0.31)** 
0.32** (0.18)** 
0.08** (0.04)** 
o** (0) ** 


ooooojo oc o 


*Sinusoidal with amplitude equal to 1.85 degrees, frequency equal to 0.16 Hz. 
** Varies with time (Figure 7); value shown is at time of interest. 


TABLE 5 - REPEAT RUNS FOR * FOR STEADY-AHEAD OPERATION 


Date of Run: ——— 18 Dec 76 21 Dec 76 ——— ———— 28 Dec 76 ——-—-_ | 30 Dec. 76 | 
Time of Run: 1100 «©1300 «=-1410 | 1740 =:1052,—s«é1613': |:«11200Ss«1310—Ss«1315~—S «1335 | 1000 | 
Run No.: 45 50 55 57 65 92 114 115° 116 117 128 
seat aan 14.08 1408 1408 | 1408 1408 14.01 | 14.08 1408 14.02 1408 | 14.00 
bir v 11.04 11.04 11.03! 11.07 11.05 11.04 | 1097 1097 1097 1097 | 10.97 
n Amplitudes (F ,)** (Ib) 
0 4380 4392 4408 | 4520 4635 4528 | 4321 4.495 4.434 4.457 | 4.528 
1 1.628 1.631 1.632 | 1532 1.631 1.622 | 1.563 1.608 1.594 1.588 | 1.610 
2 0.306 0.303 0.305 | 0340 0.335 0.330 | 0.255 0.257 0.261 0.255 | 0.292 
3 0.215 0.210 90.211 | 0.221 0.228 0.234 | 0.201 0.203 0.204 = 0.201 '||_:«0.229 
4 0.115) 0.117) 6.117 | 0.110 0.118 0.123 | 0.119 0.123. 0.123 0.121 | (0.127 
5 0.060 0.061 0.059 | 0.054 0.063 0.069 | 0.054 0.055 0.051 0.052 | 0.056 
6 0.028 0.027 0.024 | 0.028 0.031 0.033 | 0.021 0.025 0.022 0.021 | 0.027 | 
7 0.008 0.010 0.008 | 0.007 0.003 0.002 | 0.003 0.002 0.002 0.003 | 0.005 
8 0.015 0.016 0.017 | 0.008 0.010 0.008 | 0.010 0.010 06.011 0.011 | 0.013 | 
9 0.009 0.009 0.011 | 0.010 0.002 0.004 | 0.003 0.003 0.004 0.005 | 0.008 | 
10 0.008 0.007 0.009 | 0.009 0.010 0.011 | 0.010 0.009 0.009 0.010 | 0.010 
" 0.008 0.009 0.009 | 0.012 0.011 0.011 | 0012 0.014 0.014 0.016 | 0.012 
12 0.013 0.013 0.014 | 0.015 0.014 0.015 | 0.017 0.016 0.017 0.019 | 0.018 | 
13 0.009 0.009 0.010 | 0.012 0010 0.011 | 0.012 0.011 0.011 0.014 | 0.014 
14 0.006 0.006 0.007 | 0.008 0.007 0.007 | 0.008 0.006 0.005 0.008 | 0.008 | 
15 0.005 0.004 0.006 | 0.003 0.002 0.004 | 0.005 0.006 0.005 0.005 | 0.006 | 
16 0.007 0.006 0.008 | 0.005 0.003 0.003 | 0.005 0.005 0.006 0.006 | 0.006 


Phases (b,)** (deg) 


1 188 8 86118.6 119.0 | 115.3 116.8 116.5 116.7. 116.5 116.3 116.2 118.7 
2 12.0 11.2 11.6 3.5 -0.0 -0.2 1.2 1.3 3.4 3.2 -2.7 
3 38.2 37.3 37.9 19.1 24.1 23.8 21.6 21.4 21.7 20.7 26.0 
4 37.7 36.8 36.6 17.6 20.9 19.2 24.8 23.6 26.5 23.7 30.3 
5 31.9 30.4 30.4 -0.6 10.5 9.5 13.5 12.6 15.5 12.2 18.7 
6 49.7 46.5 39.8 46 13.5 16.4 11.2 5.8 6.2 3.0 35.5 
7 167.9 -169.0 -166.5 | 107.7 113.7. 116.5 | -159.1 -92.1 1610 -1688 101.7 
8 -126.4 -125.5 -121.0 |-1716 -155.3 179.3 | -137.8 -123.6 -127.9 -140.6 | -147.2 
9 -135.2 -126.3 -133.4 140.1 171.200 136.1 | -157.9 -1139 -1208 127.9 172.0 
10 172.9 -1783 -176.1 107.8 128.3 125.0 130.7. 123.6 119.1 120.3 143.5 
1 25.3 6161.8 176.2 | 106.7 1306 111.9 127.5 118.4 122.3 121.3 146.1 
12 141.7. 140.9 143.8 91.5 112.1 100.1 114.5 115.0 114.3 122.1 143.7 
13 133.7 135.6 143.2 95.4 114.5 1083 123.8 = 121.2 119.7 123.2 132.1 
14 116.9 112.7 117.6 75.5 113.2 97.0 84.3 78.3 85.8 85.0 90.6 
15 90.6 98.8 113.7 | 110.8 96.1 74.4 60.3 58.2 36.6 56.8 26.1 
16 719 78.4 88.6 55.1 746 61.0 $3.8 47.7 68.3 85.6 92.8 


. “Run 115 used for detailed analysis 
**Raw data without corrections for interactions or downstream body 
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TABLE 6 - CENTRIFUGAL AND GRAVITATIONAL LOADS* 


(Ib) (in.-Ib) (Ib) (in.-Ib) (Ib) (in.-Ib) (r/s) 

Mean, (F, M) 0.036 0.191 0.774 0.0 7.300 -0.510 14.08 
First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 
First Harmonic Phase, (Pp 14) ,, 0.0 0.0 -96.0 96.0 -159.0 0.0 

(deg) 
Mean, (F, M) 0.035 0.184 0.744 0.0 6.930 -0.485 13.76 
First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 
First Harmonic Phase, (be, um): 0.0 0.0 -96.0 96.0 -159.0 0.0 

(deg) 
Mean, (F, M) 0.030 0.161 0.642 0.0 5.860 -0.410 12.65 
First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 
First Harmonic Phase, (P¢ 14),, 0.0 0.0 -96.0 96.0 -159.0 0.0 

(deg) 
Mean, (F, M) 0.024 u.127 0.493 0.0 4.400 -0.308 10.96 
First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 
First Harmonic Phase, (¢ yy) ,. 0.0 0.0 -96.0 96.0 -159.0 0.0 

(deg) 
Mean, (F, M) 0.016 0.108 0.434 0.0 3.820 -0.207 10.21 
First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 
a Harmonic Phase, ($¢ M)4 P 0.0 0.0 ~96.0 96.0 -159.0 0.0 

deg) 3 


*The results shown here are for an Aluminum model propeller, Pp = 5.44 lbf -s?/#r* and were 
obtained by fairing the experimental air-spin data using the following restraints: 


(F, M)/n? = constant 
(F,, M4) =constant 
%,  =constant 


Except where indicated in Table &, al! other centrifugal and gravitational loads presented in this 
report are for a Bronze propellr, Pp = 14.48 Ibf-s?/tt* 
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TABLE 7 - SUMMARY OF CIRCUMFERENTIAL VARIATION OF LOADS AT THE 
SELF PROPULSION CONDITION; V = 6.52 KNOTS, 
n = 14.08 REVOLUTIONS PER SECOND 


Maximum Values Minimum Values First Harmonic Values 
(F, Max (F, Msniny (F,M), 
ee Omax dmin ean >, 
|F,M| |F, MI (F,M) 

Fea 1.43 124 0.51 308 0.42 116 
My 1.42 128 0.52 312 0.41 120 
Fy 1.40 120 0.54 292 0.38 105 
Mx, -1.34* 124 -0.61 300 0.32 294 
Foe -7.88* 252 6.30 68 6.87 54 
M2, -1.30* 204 -0.63 64 0.31 41 
Mo a 1.37 132 0.59 312 0.35 124 
Mo 44 1.36 132 0.60 312 0.34 128 
se 1.42 124 0.53 308 0.41 116 
My 1.40 128 0.54 312 0.39 120 
Fy 1.14 132 0.83 304 0.12 123 
M, -1.26* 136 -0.70 308 0.23 304 
F, 1.04 76 0.96 264 0.03 89 
M, -1.18* 204 -0.77 64 0.19 41 
Mo3 1.45 136 0.48 312 0.43 124 
Mo.4 1.60 140 0.33 312 0.56 126 


*The maximum values shown are the values which have the largest absolute value. For maximum values 
shown with negative sign, the corresponding time average values are negative using the adopted 
convention shown in Figure 1. 


TABLE 8 - TIME AVERAGE LOADS FOR STEADY~AHEAD OPERATION AT THE 
SELF PROPULSION CONDITION; V = 6.52 KNOTS, 
n = 14.08 REVOLUTIONS PER SECOND 


Total with Total with Hydrodynamic Total with Hydrodynamic 
Aluminum Bronze Loads Bronze Loads 
G'ades Blades Only Blades Only 


(1b) 3.722 3.781 3.686 é 0.0437 
(in.-Ib) 10.903 11.220 10.712 : : 0.0154 
(Ib) 3.456 4.742 2.682 ; 0.0318 
(in.-Ib)-7.765 = 7.765. = -7.765 : -0.0112 
(Ib) 7.145  —-19.279 -0.154 : -0.0018 
(in.-Ib)  -2.608 = -3.456 -2.098 : . -0.0030 
03 (in.-Ib) «6.659 5.563 7.322 eB 0.0106 


Mo4 (in.-Ib) 4.653 3.314 5.609 A 0.0081 
Mos (in.-Ib) 2.790 1.079 ie 0.0055 


For Aluminum blades pp, = 5.44 Ibf-37/#x4 
For Bronze blades pp = 14.48 Ibf-s”/x4* 


Kp = F/(pn?D*) and Ky = M/(pn2D®) are based on the density of water 
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